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In 1979 Louis Witten demonstrated that stationary axially symmetric Ein- 
stein field equations and those for static axially symmetric self-dual SU(2) gauge 
fields can both be reduced to the same (Ernst) equation. In this paper we 
use this result as point of departure to prove the existence of the mass gap 
for quantum source-free Yang-Mills (Y-M) fields. The proof is facilitated by 
results of our recently published paper, JGP 59 (2009) 600-619. Since both 
pure gravity, the Einstein-Maxwell and pure Y-M fields are described for axi- 
ally symmetric configurations by the Ernst equation classically, their quantum 
descriptions are likely to be interrelated. Correctness of this conjecture is suc- 
cessfully checked by reproducing (by different methods) results of Korotkin and 
Nicolai, Nucl.Phys.B475 (1996) 397-439, on dimensionally reduced quantum 
gravity. Consequently, numerous new results supporting the Faddeev-Skyrme 
(F-S) -type models are obtained. We found that the F-S-like model is best 
suited for description of electroweak interactions while strong interactions re- 
quire extension of Witten's results to the SU(3) gauge group. Such an extension 
is nontrivial. It is linked with the symmetry group SU(3)xSU(2)xU(l) of the 
Standard Model. This result is quite rigid and should be taken into account in 
development of all grand unified theories. Also, the alternative (to the F-S-like) 
model emerges as by-product of such an extension. Both models are related to 
each other via known symmetry transformation. Both models possess gap in 
their excitation spectrum and are capable of producing knotted/linked config- 
urations of gauge/gravity fields. In addition, the paper discusses relevance of 
the obtained results to heterotic strings and to scattering processes involving 
topology change. It ends with discussion about usefulness of this information 
for searches of Higgs boson. 
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1 Introduction 



1.1 General remcirks 

History of physics is full of situations when experimental observations lead to 

deep mathematical results. Discovery of Yang-Mills (Y-M) fields in 1954 [1] falls 
out of this trend. Furthermore, if one believes that theory of these fields makes 
sense, they should never be directly observed. To make sure that these fields 
do exist, it is necessary to resort to all kinds of indirect methods to probe them. 
Physically, the rationale for the Y-M fields is explained already in the original 
Yang and Mills paper [1]. Mathematically, such a field is easy to understand. 
It is a non Abelian extension of Maxwell's theory of clectromagnctism. In 
1956 Utiyama [2] demonstrated that gravity, Y-M and electromagnetism can 
be obtained from general principle of local gauge invariance of the underlying 
Lagrangian. The explicit form of the Lagrangian is fixed then by assumptions 
about its symmetry. For instance, by requiring invariance of such a Lagrangian 
with respect to the Abelian U(l) group, the functional for the Maxwell field 
is obtained, while doing the same operations but using the Lorentz group the 
Einstein-Hilbert functional for gravitational field is recovered. By employing 
the SU(2) non Abelian gauge group the original Y-M result [1] is recovered. 

Only Maxwell's electromagnetic field is reasonably well understood both at 
the classical and quantum level. Due to their nonlinearity, the Y-M fields are 
much harder to study even at the semi/classical level. In particular, no classical 
solutions e.g. solitons (or lumps) with finite action are known in Minkowski 
space-time. This result was proven by many authors, e.g. see [3-4] and rcfr- 
erences therein. The situation changes dramatically in Euclidean space where 
the self-duality constraint allows to obtain meaningful classical solutions [5,6]. 
These are helpful for development of the theory of quantum Y-M fields. Such 
solutions are useful in the fields other than quantum chromodynamics (QCD) 
since the self-duality equations are believed to be at the heart of all exactly 
integrable systems [7] . Although the self-duality equations originate from study 
of the Y-M functional, not all solutions [6] of these equations are relevant to 
QCD. In this paper we discuss the rationale behind the selection procedure. 
In QCD solutions of self-duality equations, known as instantons, are describ- 
ing tunneling between different QCD vacua [8]. It should be noted though 
that treatment of instantons in mathematics [9-11] and physics literature [8] 
is different. This fact is important. It is important since one of the major 
tasks of nonperturbative QCD lies in developing mathematically correct and 
physically meanigful description of these vacua. According to a point of view 
existing in physics literature the QCD has a countable infinity of topologically 
different vacua. Supposedly, the Faddeev-Skyrme (F-S) model is designed for 
description of these vacua. If this model can be used for this purpose, then 
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each vacuum state is expected to be associated with a particular knot (or hnk) 
configuration. Under these conditions the instantons are beheved to be weU 
localized objects interpolating between different knotted/linked vacuum config- 
urations [12-16]. These configurations upon quantization are expected to posses 
a tower of excited states. Whether or not such a tower has a gap in its spec- 
trum or the spectrum is gapless is the essence of the millennium prize problerrf^. 
Originally, the above results were obtained and discussed only for SU(2) gauge 
fields [17]. They were extended to SU(N) case, > 2, only quite recently [180. 
Although such a description of QCD vacua is in accord with general principles 
of instanton calculations [8] , it is in formal disagreement with results known in 
mathematics [9-11]. Indeed, it is well known that complement of a particular 
knot in is 3- manifold. Since instantons "live" in R"' (or any Riemannian 
4-manifold allowing an anti self -dual decomposition of the Y-M field (e.g. see 
Ref.[9], pages 38-39(1), this means that all knots in R^(or S"*) are trivial and 
one should talk about knotted spheres instead of knotted rings [20] . This known 
topological fact is in apparent contradiction with results of [13-15]. In this work 
we shall provide evidence that such a contradiction is only apparent and that, 
indeed, knotted configurations in are consistent with the notion of instan- 
tons as formulated in mathematics. This is achieved by using results by Floer 
[21]. It should be noted, though, that known to us "proofs" [22-24] of the 
existence of the mass gap in pure Y-M theory done at the physical level of rigor 
ignore instanton effects altogether. Among these papers only Ref. [22] uses the 
F-S SU(2) model for such mass gap calculations. It also should be noted that 
results of such calculation sensitively depend upon the way the F-S model is 
quantized. For instance, in the work by Faddeev and Niemi [25], done for the 
SU(2) gauge group, the results of quantization produce gapless spectrum. To fix 
the problem the same authors suggested to extend the original model in ad hock 
fashion. Other authors, e.g. see Ref. [26], proposed different ad hoc solution of 
the same problem. 

The above results are formally destroyed by the effects of gravity. Indeed, in 
1988 Bartnik and McKinnon numerically demonstrated [27] that the combined 
Y-M and gravity fields lead to a stable particle-like (solitonic) solutions while 
neither source- free gravity nor pure Y-M fields are capable of producing such so- 
lutionfl Such situation has interesting cosmological ramification^ [28] causing 
disappearance of singularities in spacetime as shown by Smoller et al [29] . In 
this work we do not discuss implications of these remarkable results. Instead, in 
the spirit of Floer's ideas [21], we argue that even without taking these results 
into account, the effects of gravity on processes of high energy physics are quite 
substantial. 

^E.g. see 

http:/ /www. claymath.org/millennium/Yang-Mills_Theory/' 

■^In this work, in accord with experimental evidence, we demonstrate that N < 3. 
^In physics literature, both anti and self dual instantons are allowed to exist, e.g. see 
Ref. [19], page 481. 

^More accurately, neither pure Y-M fields nor pure gravity have nontrivial static globally 
regular (i.e.nonsingular, asymptotically fiat) solitons. 
^E.g. Einstein- Y-M hairy black holes 
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1.2 Statements of the problems to be solved 

In this paper several problems arc posed and solved. In particular, we would like 
to investigate the physics and mathematics behind gravity- Y-M correspondence 
discovered by Louis Witten [30] for SU(2) gauge fields. Is this correspondence 
accidental? If it is not accidental, how it should be related to commonly shared 
opinion that the Standard Model (SM) of particle physics does not account 
for gravity? Can this correspondence be extended to other gauge fields, e.g. 
SU(N), N>2 ? If the answer is "yes", will such correspondence be valid for all 
N's or just for few? In the last case, what such a restriction means physically? 
How the noticed correspondence is helping to solve the gap problem? What role 
the F-S model is playing in this solution? Is this model instrumental in solving 
the gap problem or are there other aspects of this problem which the F-S model 
is unable to account? How this correspondence affects known string-theoretic 
and loop quantum gravity (LQG) results ? What place the topology-changing 
(scattering) processes occupy in this correspondence? Is there any relevance of 
the results of this work to searches for Higgs boson? 

1.3 Organization of the rest of the paper and summeiry of 
obtained results 

Sections 2,3 and 6, and Appendix A are devoted to detailed investigation of 
gravity- Y-M correspondence. Section 4 is devoted to the physics-style exposition 
of works by Andreas Floor [11, 21] on Y-M theory with purpose of connecting his 
mathematical formalism for Y-M fields with the F-S model. In the same section 
we also consider the Y-M fields monopole and instanton solutions and their 
meaning and place within Floor's theory. Om exposition is based on results of 
Sections 2 and 3. Section 5 is entirely devoted to solution of the gap problem for 
pure Y-M fields. Although the solution depends on results of previous sections, 
numerous additional facts from statistical mechanics and nuclear physics are 
being used. In Section 6 we discuss various implications/corollaries of the 
obtained results, especially for the SM of particle physics. In Section 7 we 
discuss possible directions for further research based on the results presented 
in this paper. These include (but not limited to): connections with the LQG, 
the role and place of the Higgs boson, relationship between real space-time 
scattering processes of high energy physics and processes of topology change 
associated with such scattering. Based on the results of this paper, we argue 
that this task can be accomplished with help of the formalism developed by G. 
Perelman for his proof of the Poincare' and geometrization conjectures. 
The major new results of this paper are summarized as follows. 

1. In subsection 5.4.4, while solving the gap problem, we reproduced by 
employing entirely different methods, the main results of the paper by Korotkin 
and Nicolai [31] on quantizing dimensionally reduced gravity. From these results 
it follows that for gravity and Y-M fields possessing the same symmetry the 
nonperrturbative quantization proceeds essentially in the same way. 

2. In subsection 6.3 we demonstrated that gravity- Y-M correspondence dis- 
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covered by L.Wittcn for gauge group SU(2) can be extended only to the SU(3) 
gauge group. This group contains SU(2)xU(l) group as a subgroup. This fact 
allowed us to come up with the anticipated (but never proven!) conclusion about 
symmetry of the SM. It is given by SU(3)xSU(2)xU(l). The obtained result is 
very rigid. It is deeply rooted into not widely known/appreciated (discussed in 
Appendix A) properties of the gravitational field. It is these properties which 
ultimately determine the conditions of gravity- Y-M correspondence. 

3. The latest papers Refs. [32-34] are aimed at reproduction of the classifi- 
cation scheme of particles and fields in the SM within the framework of LQG 
formalism. These results match perfectly with the results of our paper be- 
cause of the noticed and developed gravity- Y-M correspondence. In view of this 
correspondence, the results of Refs. [32-34] can be reproduced with help of min- 
ima] gravity model described in subsections 3.2, 3.4, and 7.2 . This minimal 
model has differential-geometric /topological meaning in terms of the dynamics 
of the extended Ricci fiow [35,36]. Such a fiow is the minimal extension of the 
Ricci flow now famous because of its relevance in proving the Poincare' and 
geometrization conjectures. 

4. The formalism developed in this paper explains why using pure gravity 
one can talk about the particle/field content of the SM. Not only it is compatible 
with just mentioned LQG results but also with those, coming from noncommu- 
tative geometry [37], where it is demonstrated that use of pure gravity (that is 
"minimal model") combined with 0- dimensional internal space is sufficient for 
description of the SM. 

2 Emergence of the Ernst equation in pure grav- 
ity and Y-M fields 

2.1 Some facts about the Ernst equation 

Study of static vacuum Einstein fields was initiated by Weyl in 1917. Con- 
siderable progress made in later years is documented in Ref.[38]. To develop 
formalism of this paper we need to discuss some facts about these static fields. 
Following Wald [39] , a spacetime is considered to be stationary if there is a one- 
parameter group of isometrics at whose orbits are time-like curves e.g. see [40]. 
With such group of isometrics is associated a time-like Killing vector Fur- 
thermore, a spacetime is axisymmetric if there exists a one-parameter group of 
isometries x,p whose orbits are closed spacelike curves. Thus, a spacelike Killing 
vector field tp^ has integral curves which are closed. The spacetime is station- 
ary and axisymmetric if it possesses both of these symmetries, provided that 
o-* °X<j, = X<j> oo-t • If C = im) ^ = (^) that [^, -0] = 0, one can choose 
coordinates as follows: =t,x^ = (p, x"^ = p,x^ = z. Under such identification, 
the metric tensor g/j^^, becomes a function of only a;^and x^. Explicitly, 

ds'^ = -V{dt - wd(f>f + V-\p'^d(l? + e2T(rfp2 + dz'^% (2.1) 
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where functions V, w and 7 depend on p and z only. In the case when V = 
l,w — 'J — 0, the metric can be presented as ds^ = — (dt)^ +{ds)'^, where 

{dSf = p^dcj)^ + dp^ + dz^ (2.2) 

is the standard flat 3 dimensional metric written in cylindrical coordinates. The 
four-dimensional set of vacuum Einstein equations Rij — with help of metric 
given by Eq.(2.2) acquires the following form 

V -{V-^VV + p-^V^wVw} (2.3a) 

and 

V • {p-^V^Vw} ^ 0. (2.3b) 

In these equations V- and V are three-dimensional flat (that is with metric given 
by Eq.(2.2)) divergence and gradient operators respectively. In addition to these 
two equations, there are another two needed for determination of factor 7 in the 
metric, Eq.(2.1). They require knowledge of V and w as an input. Solutions 
of Eq.s(2.3) is described in great detail in the paper by Reina and Trevers [41] 
with final result: 

(Ree)V^e = Ve- Ve. (2.4) 

This equation is known in literature as the Ernst equation. The complex po- 
tential e is defined mhy e — V + iu with V defined as above and lo being an 
auxiliary potential whose explicit form we do not need in this work. As it was 
recognized by Ernst [42,43] such an equation can be also used for description of 
the combined Einstein-Maxwell fields. We shall exploit this fact in Section 6. 
In Appendix A and in Section 6 we provide proofs that knowledge of static vac- 
uum solutions of the Ernst equation is necessary and sufficient for restoration of 
static Einstein-Maxwell fields!^ Fields other than Y-M should be also restorable 
0. To proceed, we need to list several properties of the Ernst equation to be 
used below. First, following [41] and using prolate spheroidal coordinates, the 
Ernst equation reproduces the Schwarzschild metric, and with another choice 
of coordinates it reproduces the Kerr and Taub-NUT metric. Thus, the Ernst 
equation is the most general equation describing physically interesting vac- 
uum spacetimes compatible with the Cauchy formulation of general relativity 
[39,40,44,45]. Such a formulation is convenient staring point for quantization 
of gravitational field via superspace formalism [39] leading to the Wheeler- De 
Witt equation, etc. Since in ttiis work we advocate different approacli to quan- 
tization of gravity, this topic is not being discussed further. Second, following 
Ref.[38], page 283, a stationary solution of Einstein's field equations is called 
static if the timelike Killing vector is orthogonal to the Cauchy surface. In such 
a case from the Table 18.1. of the same reference it follows that the Ernst 

^Surprisingly, upon changes of variables in these static solutions, the exact results for 
propagating gravitational waves can be obtained as well. 

^This is so because each of these fields is a source of gravitational field which, in turn, can 
be eliminated locally. See Appendix A. 
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potential e is real. This observation allows us to simplify Eq.(2.4) considerably. 
For the sake of notational comparison with Ref.[38] we redefine the potential 
e = V + In the static case we have e = —F = —e^^. Using this result in 
Eq.(2.4) produces 

^p,zu = 0, (2.5) 

where Ap^^ is flat Laplacian written in cylindrical coordinates defined by the 
metric, Eq.(2.2). 

2.2 Isomorphism between the SU(2) self-dual gauge and 
vacuum Einstein field equations 

This isomorphism was discovered by Louis Witten in 1979 [30]. His work was 
inspired by earlier works of Ernst [42] and Yang [46]. To our knowledge, since 
time when Ref. [30] was published such an isomorphism was left undeveloped. 
In this paper we correct this omission in order to demonstrate that when both 
fields are mathematically indistinguishable, their quantization should proceed 
in the same way. The result analogous to that discovered by Witten was ob- 
tained using different arguments a year later by Forgacs, Horvath and Palla [47] 
and, in a simpler form, by Singleton [48]. These authors used essentially the 
paper by Manton, Ref. [49], in which it was cleverly demonstrated that the 't 
Hooft-Polyakov monopole can be obtained without actual use of the auxiliary 
Higgs field. Both Refs. [47,48] and the original paper by Witten [30] use the 
axial symmetry of either gravitational or Y-M fields essentially. Only in this 
case it can be shown that the axisymmetric version of the self-duality equations 
obtained by Manton can be rewritten in the form of the Ernst equation. In 
the light of above information, following Ref. [5], we shall discuss briefly con- 
tributions of Yang and Witten. For this purpose, we need to consider first the 
following auxiliary system of linear equations 

= X*;*t = T*. (2.6) 

Here "if^ = ^'^d — In. this system X and T are square matrices 

of the same dimension and such that 

Xt-T, + [X,T] =0 (2.7) 

This result easily follows from the compatibility condition: 'ifxt = "^tx- The 
matrices X and T can be realized as 

(.-,:<) V)'-<c 'a) <-) 

with ^ being a spectral parameter and. A, B and C being some Laurent poly- 
nomials in C. The above system can be extended to four variables xi, a;2, ii, ^2 



°The minus sign in front of F is written in accord with conventions of Chapter 30.2 of the 
1st edition of Ref. [38] . 
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in a simple minded fashion as follows 



(|_+C^)* = (T,+iT.)*. (2.9b) 

In the most general case, the matrices Xi,X2,Ti,T2 are made of functions 
which "live" in C*. They are representatives of the Lie algebra sl{n, C) of n x n 
trace-free matrices. The compatibility conditions for this case arc cqnivalcnt to 
the self-duality condition for the Y-M fields associated with algebra sl{n, C).It 
is instructive to illustrate these general statements explicitly. 

In the (anti)sclf-duality condition for the Y-M curvature reads: *F 
i so that for the self-dual case we obtain: 

Foi = F23, F02 = F3i,Fo3 = F12. (2.10) 

In the "light cone" coordinates cr = '^i^i + 1x2), t = -^{^a + '1x3) the Y-M 
field one- form can be written as A/j^dx^ = A^da + A^dr + Agda + Afdf with the 
overbar labeling the complex conjugation. In such notations Aq = -^{A^ + Af), 

Ai = -^{A^+A^), A2 = -^{A„ - A^), A3 = -^{Ar - Af). In these notations 
Eq.s (2.9) acquire the following form 

F^r = 0, F^f = and F^^ + Frf = 0. (2.11) 

They can be obtained as compatibility condition for the isospectral linear prob- 
lem 

[d, + Qdr)^ = {A„HM)'if and (9, - C^^)* = {A^ - C^^)*, (2.12) 

where the spectral parameter is C, and * is the local section of the Y-M fiber 
bundle. The compatibility condition reads: {d^ — (df){da + C^r)* = {da + 
(df){da — (df)^, thus leading to 

[Far - C{Faa + Frf) + C'f'af ]* = 0. (2.13) 

This equation allows us to recover Eq.s(2.11). The first two equations of 
Eq.s(2.11) can be used in order to represent the A -fields as follows: Aa- = 
{daC)C-\ Ar ^{drC)C-\ Aa={dgD)D-^ and Af = {drD)D-^, where 
both C and D are some matrices in the Lie group G, e.g. G = SU{2). By 
introducing the matrix M = C~^D G G the last of equations in Eq.(2.11) 
becomes 

da{M-'^daM) + df{M-^drM) = 0. (2.14a) 

Thus, the self-duality conditions for the Y-M fields are equivalent to Eq.(2.14a). 

For the future use, following Yang [46], we notice that in such formalism the 
gauge transformations for Y-M fields are expressible through D — > DE and 
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C ^ CE so that F^g. E-'^F^^E and F^r E-'^FrfE with the matrix 
E = E{a,a,T,f) S SU{2) leaving self-duality Eq.s(2.10) (or (2.13)) unchanged. 

To connect Eq. (2.14a) with the Ernst equation, following L.Wittcn [30] it 
is sufficient to assume that the matrix M is a function of p = \/xl + x\ and 
2; = a;3. In such a case it is useful to remember that = 2a a and z = :^(''"~^)- 
With help of these facts Eq.(2.14a) can be rewritten as 

dpipM-'^dpM) + pd^{M-^d^M) = 0. (2.14b) 

By assuming that the matrix M is representable by the S'L(2, i?)-type matrix, 
and writing it in the form 

Eq.(2.14b) is reduced to the pair of equations 

VV^V = Vy • - V$ • V$ and yV^$ = 2VV • v$. 

With help of the Ernst potential e = V + i^ these two equations can be brought 
to the canonical form of the Ernst equation, Eq.(2.4). Below, we shall provide 
sufficient evidence that such a reduction of the Ernst equation is compatible with 
analogous reduction in instanton/monopole calculations for the Y-M fields. 

3 From analysis to synthesis 
3.1 General remEirks 

The results of previous section demonstrate that for axially symmetric fields 
both pure gravity and pure self-dual Y-M fields are described by the same 
(Ernst) equation. In this section we reformulate these results in terms of the 
nonlinear sigma model with purpose of using such a reformulation later in the 
text. To do so we need to recall some results from our recent works, Ref.s[50,51]. 
In particular, we notice that under conformal transformations g = e^'^g in d- 
dimensions the curvature scalar R{g) changes as follows: 

R{g) = e-^^{R{g) - 2{d - l)AgU - (rf - l){d - 2) \Vguf}- (3,1) 

Since this equation is Eq.(2.11) of our Ref.[50] we shall be interested only in 
transformations for which R{g) is a constant. This is possible only if the total 
volume of the system is conserved. Under this constraint we need to consider 
Eq.(3.1) for d = 3 in more detail. Without loss of generality we can assume that 
initially R{g) = 0. For this case we shall write g = go so that Eq.(3.1) acquires 
the form 

Rig) = -2e-2"[2A,„w + IVgo^f] (3-2) 

in which A^,, is the flat space Laplacian. Now we can formally identify it 
with that in Eq.(2.5). Accordingly, we shall be interested in such conformal 
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transformations for which Ag„u = in Eq.(3.2). If they exist, Eq.(3.2) can be 
rewritten as 

e'-m = -2 (v,„^^) • (Vgo") • (3.3) 
This allows us to interpret Eq.(3.3) and 

Ago" = (3.4) 

as interdependent equations: solutions of Eq.(3.4) determine the scalar cur- 
vature R{g) in Eq.(3.3). Clearly, under conditions at which these results are 
obtained only those solutions of Eq.(3.4) should be used which yield the con- 
stant scalar curvature R{g). Eq.(3.3) contains information about the Ricci 
tensor. To recover this information we notice that gij — —e^^Sij. Therefore we 
obtain: 

R,j{g)^2V,uVjU, (3.5) 

in accord with Eq.(18.55) of Ref.[38] where this result was obtained by employing 
entirely different arguments. From the same reference we find that Eq.(3.4) 
comes as result of use of the contracted Bianci identities applied to i?ij(gjj. 
It is instructive to place the obtained results into broader context. This is 
accomplished in the next subsection. 



3.2 Connection with the nonlinear sigma model 

Some time ago Neugebauer and Kramer (N-K), Ref.[38], obtained Eq.s(3.4) and 
(3.5) using variational principle. In less general form this principle was used pre- 
viously by Ernst [42] resulting in now famous Ernst equation. Neugebauer and 
Kramer proposed the Lagrangian and the associated with it action functional 
Sn-k producing upon minimization both Eq.s(3.4) and (3.5). To describe these 
results, we also use some results by Gal'tsov [52]. 
The functional Sn-k is given by 

Sn-k = \j V§[R{9) ~ 9''GAB{^)d,v^d,^'']dPx, (3.6) 

M 

easily recognizable as three-dimensional nonlinear sigma model coupled to 3-d 
Euclidean gravity. The number of components for the auxiliary field ip as well 
as the metric tensor GABi'p) of the target space is determined by the problem 
in question. In our case upon variation of Sn-k with respect to ipf and gij we 
should be able to re obtain Eq.s(3.4) and (3.5). To do so, following Ref. [53], we 
introduce the current 

J, = AT^diM. (3.7) 

In view of results of subsection 2.2, we have to identify the matrix M with that 
defined by Eq.(2.15) and, taking into account Eq.(2.14a), the index i should 

9Ref.[38], page 283, bottom 
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take two values: tr and r. With such definitions we can replace the functional 

Sn-k by 

s = ^Tmig) - a'^-^HJi Jk)]d'x. (3.8) 

M 

The actual calculations with such type of functionals can be made using 
results of Ref.[53]. Thus, using this reference we obtain, 

Rijig) = \tr{Ji Jj) (3.9) 

and 

diJi =0. (3.10) 

Evidently, by construction Eq.(3.10) coincides with Eq.(2.14a) and, ultimately, 
with Eq.(3.4). It is also easy also to check that Eq.(3.9) does coincide with 
Eq.(3.5). For this purpose it is sufficient to notice that 

tr{Ji Jj) = -tr{diMdjM-'^). (3.11) 

To check correctness of our calculations the entries of the matrix M, Eq.(2.15), 
can be restricted to V (that is we can put $ = 0). Since V = —F = — e^" (e.g. 
see discussion prior to Eq.(2.5)), a simple calculation indeed brings Eq.(3.9) 
back to Eq.(3.5) as required. 

It is interesting and important to observe at this point that the equa- 
tion of motion, Eq.(3.10), formally is not affected by effects of gravity. This 
conclusion requires some explanation. From subsection 2.2, especially from 
Eq.s(2.14),(2.15), it should be clear that Eq.(3.10) is the Ernst equation deter- 
mining gravitational field. Hence, it is physically wrong to expect that it is 
going to be affected by the effects of gravity. Eq.s(3.9) and (3.10) are the same 
as Eq.s(3.5) and (3.4) whose meaning was explained in the previous subsection. 
Clearly, the functional, Eq.(3.8), can be used for coupling of other fields to 
gravity. This is indeed demonstrated in Ref.[52]. This is done with purpose 
of connecting results for the nonlinear sigma models with those for heterotic 
strings. We would like to discuss this connection now since it will be used later 
in the text. 

3.3 Connection with heterotic string models 

The functional <S, Eq.(3.8), is related to that for the heterotic string model, 
e.g. see Ref . [54] . For such a model the sigma model-like functional is obtainable 
from 10 dimensional supersymmetric string model by means of compactifica- 
tion scheme (ideologically similar to that used in the Kaluza-Klein theory of 
gravity and electromagnetism) aimed at bringing the effective dimensionality 
to physically acceptable values (e.g. 2, 3 or 4). For dimensionality D <10 such 
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compactified/reduced action functional reads (e.g. see Ref.[54], Eq.(9.1.8)): 



-i(M-i),,F;,F^>'- + hr{d^Md^^M-% (3.12) 

The compactification procedure is by no means unique. There are many ways 
to make a compactified action to look exactly like that given by Eq.(3.8) (e.g. 
see [55]). Evidently, there should be a way to relate such actions to each other 
since they all are having the same origin - 10 dimensional heterotic string action. 
Because of this, we would like to make some comments on action given by 
Eq.(3.12) by specializing to I? = 3 for reasons explained in Refs[68,69] and 
to be clarified below, in Section 6. Under such conditions if we require the 
dilaton 0, the antisymmetric iJ-field (associated with string orientation) and 
the electromagnetic field F to vanish, the remaining action will coincide with 
that given by Eq.(3.8). Because of this, the following steps can be made. 

First, as explained in our work, Ref. [50] , for closed 3-manifolds we can/should 
drop the dilaton field (p. Second, by properly selecting string model we can ignore 
the antisymmetric field H. Third, by taking into account results of Appendix A 
we can also drop the electromagnetic field since it can be always restored from 
pure gravity. Thus, we end up with the action functional 5, Eq.(3.8), which 
we shall call "minimar. In Section 6 we shall provide evidence that its mini- 
mality is deeply rooted into gravity- Y-M correspondence which does not leave 
much room for "improvements" abundant in physics literature. We shall begin 
explaining this fact immediately below and will end our arguments in Section 
6. 



heterotic 



3.4 The extended Ricci flow 

Thus far use of the variational principle apparently had not brought us any 
new results (at least at the classical level). Situation changes in the light of 
recent work by List [35]. Following Ref.s[35,36] , it is convenient to introduce 
Perelman-like entropy functional T{gij,u, f) 

Hmi,uJ) = j{R(g)-2\Vgu\^ + \Vgf\^)e-Uv (3.13) 

M 

coinciding with Eq. (7.22b) of our work, Ref. [50], when u = 00. Because of 
this observation, if formally we make a replacement 7?.(g; u) — R{g) — 2 |Vu| 
in Eq.(3.13), we are able to identify Eq.(3.13) with Perelman's entropy func- 
tional enabling us to follow the same steps as were made in Perelman's papers 
aimed at proof of the geometrization and Poincare' conjectures. Such a program 

^"it should be noted that there is an obvious typographical error in Eq. (7.22b): the term 
|Vh/|^ is typed as |Vh/| . 
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was indeed completed in the PhD thesis by List [36] . Minimization of entropy 
functional J^{gij , /) produces the following set of equations 

d 

—gij = -2{Rij + ViVjf) +4ViuVjU, (3.14a) 
d 

-u = AgU-{Vu)-{Wf), (3.14b) 

and 

^/ = -^-Ag/ + 2|V«|% (3.14c) 

coinciding with Eq.s(7.28a), (7.28b) of our work, Ref.[50], when u = O.In these 
equations IVgw]^ = g^^ViuVjU , etc. From the next section and results below it 
follows that physically we should be interested in closed 3 manifolds. For such 
manifolds one can use Lemma 2.13, proven by List [36], which can be formulated 
as follows: 

Let g,u, / be a solution of Eq.s(3.14) for t G [0,T) on a closed manifold M. 
Then the evolution of the entropy is given by 

dtT{gij,u,f) = 2 J [\^Zi,{g■,u) + V^VJf\^ + 2{AgU-{Vu) ■ {Vf)f]e~Uv > 0. 

M 

(3.15) 

Thus, the entropy is non decreasing with equality taking place if and only if the 
solution of Eq.(3.14) is a gradient soliton. This happens when the following two 
conditions hold 

TZijig; u) + ViVjf = and AgU - (Vu) • (V/) = 0. (3.16) 

For M = the result of Perelman, Eq.(7.30) of Ref [50], for steady gradient 
soliton is reobtained, as required. Since for closed compact manifolds / = const 
Eq.s(3.16) coincide with Eq.s(3.4) and (3.5) as anticipated. Thus, existence 
of steady gradient solitons in the present context is equivalent to existence of 
sohitions of static Einstein's equations for pure gravity. This fact alone could 
be mathematically interesting but requires some reinforcement to be of interest 
physically. We initiate this reinforcement process in the following subsection. 

3.5 Relationship between the F-S and the Ernst function- 
als 

The F-S functional was mentioned in the Itroduction. In this subsection we 
would like to initiate study of its connection with the Ernst functional. We 
begin with the following observation. In steps leading to Eq.(2.14b) (or (3.10)) 
the Euclidean time coordinate xo was eventually dropped implying that solu- 
tions of selfduality for Y-M equations, when substituted back into Y-M action 
functional, will produce physically meaningless (divergent) results. While in 
subsection 4.4 we discuss a variety of means for removing of such apparent 
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divergence, in this subsection we notice that aheady Ernst [42] suggested the 
action functional whose minimization produces the Ernst equation. He gave two 
equivalent forms for such a functional, now bearing his name. These are either 

f Ve • Ve* 
SeA^]= dv-—^ (3.17) 
J (Ree 



M 

or 



SEAi] = I dv^r^- (3-18) 

Minimization of -S'Bje] leads to Eq.(2.4) while functional, Eq.(3.18), is obtained 
from SeJc] by means of substitution: e = (^ — l)/(^ + 1). In both functional 
dv is 3-dimensional Euclidean vohime element so that apparently the manifold 
M is just (or, with one point compactification, it is 5"^). Evidently, both 
S'eJe] and /S^aK] are functionals for the nonlinear sigma model. If we drop 
the curvature term in Eq.(3.6) such truncated functional can be identified, for 
example, with S'^aK]- This explains why Eq.(3.10) is formally unaffected by 
gravity. In mathematics literature the nonlinear sigma models are known as 
harmonic maps. Since Reina [56] demonstrated that the functional S'salC] 
describes the harmonic map from to H^, it is not too difficult to write 
analogous functional S'bsI^] describing the mapping from to 5^. It is given 

by 

SeA^] = I dv ^i'^^^ (3.19) 

J (^r+1) 

M 

and is part of the F-S model. If needed, both <S'£;2[^] and S'^si^] can be supple- 
mented by additional (topological) terms which in the simplest case are wind- 
ing numbers. Thus, we shall be dealing either with the truncated F-S model, 
Eq.(3.19), or with its hyperbolic analog, Eq.(3.18). The choice between these 
models is nontrivial and will discussed in detail in Section 6. To facilitate this 
discussion, we need to observe the following. In the static case, we argued, e.g. 
see Eq.(2.5), that e = —F = — e^". Substitution of this result back into S^Je] 
produces (up to a constant) the following result: 

Se^ [e] = j dvVu ■ Vu (3.20) 

M 

leading to Eq.(2.5) as anticipated. At the same time, consider the following H-E 
action functional 

Sh-e[9] = j dv^R{g), (3.21) 

M 

and take into account Eq.(3.3) and the fact that gij = —e^^Sij . Straightforward 
calculation leads us then to the result (up to a constant): 

SH-E[g] = - j dvVu ■ Vu. (3.22) 

M 
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The minus sign in front of the integral is important and will be explained be- 
low. Before doing so, we notice that the Ernst functional (in whatever form) 
is essentially equivalent to the H-E functional! Since in the original paper by 
Ernst M is (or S^), apparently, such a functional should be zero. This is 
surely not the case in general but the explanation is nontrivial. Suppose that 
minimization of the Ernst functional leads to some knotted/linked structure^"]. 
If such knots/links are hyperbolic then, by construction, complements of these 
knots/links in are H'^ modulo some discrete group. This conclusion is in 
accord with properties of the Ernst equation discovered by Reina and Trevers 
[41]. Following this reference, we introduce the complex space C x C = so 
that V z = (u,w) e the scalar product z*^z°' can be made with the metric 
T^ap — diag{l, —1}. Furthermore, the Ernst Eq.(2.4) can be rewritten with help 
of substitution e = (u — v)/{u + v) as the set of two equations 

z°'z*^W^z^ = 2z;Vz" • Vzf^. (3.23) 

Such a system of equations is invariant with respect to transformations from 
unimodular group SU(1,1) which is equivalent to SL(2,C). But SL(2,C) is the 
group of isometrics of hyperbolic space H'^ as was discussed extensively in our 
work, Ref.[57]. Thus, minimization of both the F-S and Ernst functionals should 
account for knotted/linked structures. This conclusion is strengthened in the 
next subsection. 



3.6 Relationship between the Ernst and Chern-Simons 
functionals 

Even though we need to find this relationship anticipating results of the next 
section, by doing so, some unexpected connections with previous subsection 
are also going to be revealed. For this purpose, we notice that for -u = the 
functional !F{gij,u, f) introduced earlier is just Perelman's entropy functional. 
As such, it was discussed in our work, Ref.[50]. Evidently, both F and Perel- 
man's functional can be used for study of topology of 3-manifolds. We believe, 
though, that use of Perelman's functional is more advantageous as we would 
like to explain now. For this purpose, it is convenient to introduce the Raleigh 
quotient \g via 

/dT/(4|V^|Vi?(g)^2) 
(3-24) 



M 



e.g. see Eq.(7.24) of [50], to be compared against the Yamabe quotient [p ~ 



and a = 4^) 



Yg = "^—r-. — Zn- = ( r /7T,„D ) / ^Vaia (Vg^) +R{g)^ } = —2 



We shall postpone detailed discussion of this topic till Section 6. 
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also discussed in [50]. Because of similarity of these two quotients the question 
arises: Can they be equal to each other? The affirmative answer to this question 
is obtained in Ref.[58]. It can be formulated as 

Theorem [58]. Suppose that 7 is a conformal class on M which does not 
contain metric of positive scalar cMrvature. Then 

= su\>XgV{g)i = A(M), (3.25a) 

where A(M) is Perelman's A invariant. Equivalently, 

^gV{9)i < Y^, (3.25b) 

where V{g) = J d'^xy/g is the volume. 

The equality happens when g is the Yamabe minimizcr. It is metric of 
unit volume for manifold M of constant scalar curvature (which, according to 
theorem above, should be negative so that M is hyperbolic 3-manifold). Only 
for hyperbolic 3-manifolds whose Yamabe invariant (M) = sup^ Y^ the 
gravitational Cauchy problem for source-free gravitational field is well posed 
[45,46]. For g which is Yamabe minimizer we have Sh~e[9] < Y^. This result 
can be further extended by noticing that AfSH-Eig] = CS{A) , where TV is some 
constant whose value depends upon the explicit form of the gauge field A, and 
CS{A) is the Chern-Simons invariant to be described in the next section. 

To demonstrate that J\fSH-E[g] = CS{A) it is sufficient to use some re- 
sults from works by Chern and Simons [59] and by Chern [60]. In [59] it was 
proven that: a) the Chern-Simons (C-S) functional CS{A) (to be defined in 
next section) is a conformal invariant of M (Theorem 6.3. of [59]) and, b) that 
the critical points of CS{A) correspond to 3-manifolds which are (at least lo- 
cally) conformally flat (Corollary 6.14 of [59]). Subsequently, these results were 
reobtained by Chern, Ref.[60], in much simpler and more physically sugges- 
tive way. In view of these facts, at least for Yamabe minimizers we obtain, 
CS{A) = AfSrif], where J\f is some constant (different for different gauge 
groups). That this is the case should come as not too big of a surprise since 
for Lorentzian 2-|-l gravity Witten, Ref.[61], demonstrated the equivalence of 
the Hilbert-Einstein and C-S functional without reference to results of Chern 
and Simons just cited. At the same time, the Euclidean/Hyperbolic 3d gravity 
was discussed only much more recently, for instance, in the paper by Gukov, 
Ref.[62]. To avoid duplications we refer our readers to these papers for further 
details. 

4 Floer-style nonperturbative treatment of Y- 
M fields 

4.1 Physical content of the Floer's theory 

Striking resemblance between results of nonperturbative treatment of 4-dimensional 
Y-M fields and two dimensional nonlinear sigma model at the classical level is 
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well documented in Ref.[63]. Zero curvature equations, e.g. Eq.(2.7), can be 
obtained either by using the two- dimensional nonlinear sigma model or three- 
dimensional C-S functional. As discussed in previous section, the self-duality 
condition for Y-M fields also leads (upon reduction) to zero curvature condition. 
Since the Ernst equation describing static gravitational (and electrovacuum) 
fields is obtainable both from conditions of self-duality for the Y-M field and 
from minimization of 3-dimensional nonlinear sigma model, it follows that 3-d 
gravitational nonlinear sigma model, Eq.(3.8), contains nonperturbative infor- 
mation about Y-M fields. Furthermore, in view of results of Appendix A, it 
also should contain information about the static electromagnetic fields, for the 
combined gravitational and electromagnetic waves and, with minor modifica- 
tions, for the combined gravitational, electromagnetic and neutrino fields. The 
nonperturbative treatment of Y-M fields is usually associated either with the in- 
stanton or monopole calculations. This observation leads to the conclusion that, 
at least in some cases, zero curvature equation should carry all nonperurbative 
information about Y-M fields. This point of view is advocated and developed by 
Floer [11,21] . Below,we shall discuss Floer's point of view now in the language 
used in physics literature. For the sake of illustration, it is convenient to present 
our arguments for Abelian Y-M (that is electroma gne tic) fields first. 
The action functional S in this case is given bjo 

t 

S^^ J dt J dv[E^ - B^], (4.1) 

M 

where B = V x A and E = — V(y3 — ^A , (p = Aq. It is known that, at least 
for electromagnetic waves, it is sufficient to put Aq = (temporal gauge). In 
such a case the above action can be rewritten as 

t 

S[A] ^^JdtJ dv[A^ - (V X A)^], (4.2) 

M 

where A =^A. From the condition SS/SA = wc obtain ^ = V x B. The 
definition of B guarantees the validity of the condition V • B = while from the 
definition of E we get another Maxwell equation ^ = — V x E. The question 
arises: will these results imply the remaining Maxwell's equation V • E = 
essential for correct formulation of the Cauchy problem? If such a constraint 
satisfied at i = 0, naturally, it will be satisfied for i > 0. Unfortunately, for t = 
the existence of such a constraint does not follow from the above equations and 
should be introduced as independent. This causes decomposition of the field 
A as A = A|| + Aj^. Taking into account that E = — ^A, we obtain as well 
V • (E|| +Ej^). Then, by design V • E_l — 0, while V • E|| remains to be defined 
by the initial and boundary data. Because of this, it is always possible to choose 
All = and to use only A_l for description of the field propagation [64]. Hence, 



Up to an unimportant scale factor. 
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the action functional S can be finally rewritten as 



S[Aj_] = \jdt j dv[kl - (V X A^)']. (4.3) 

M 

In such a form it can be used as action in the path integrals, e.g. see Ref. [64], 
page 152, describing free electromagnetic field. Such path integral can be eval- 
uated both in Minkowski and Eucldean spaces by the saddle point method. 
There is, however, a closely related method more suitable for our purposes. It 
is described in the monograph by Donaldson, Ref. [11]. Following this reference, 
we replace time variable t by —it in the functional >5'[Aj^] .Consider now this 
replacement in some detail. We havj^ 

T 

\jdT j dv[kl + (V X A^)'] 

M 
T 

= lJdT{Jdv[[A^ + {VxA^)f--^{A^-VxA^)]). (4.4) 

M 

Since variation of Aj^ is fixed at the ends of r integral, the last term can be 
dropped so that we are left with the condition 

^A^= (4.5) 

extremizing the Euclidean action 5'£;[Aj_]. The above results are transferable 
to the non Abelian Y-M field by continuity and complementarity. Since in 
the Abelian case fields E and B are dual to each other, by applying the curl 
operator to both sides of Eq.(4.5) (and removing the subscript _L) we obtain 
the equivalent form of self-duality equations in accord with those on page 33 of 
Ref. [6] . This calculation provides an independent check of Donaldson's method 
of computation. Since the (anti)self-duality condition in the Abelian case can 
be written as B = =fE [9]. and since E = — ^A, we conclude that Eq.(4.5) 
is the self-duality equation. This conclusion is immediately transferable to the 
non Abelian Y-M case where the analog of Eq. (4.5) is 

|-A = *F(A(t)), (4.6) 

in accord with Floer. The symbol * denotes the Hodge star operation in 3 di- 
mensions. Following Donaldson [11] this result should be understood as follows. 

3 

Introduce a connection matrix A = Aadr + ^ Aidxi such that both Aq and Ai 

i=l 

depend upon all four variables t, xi,X2 and 2:3. In the temporal gauge Aq should 



'We shall assume (without loss of generality) that Ax is coUinear with Aj 
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be discarded so that t becomes a parameter in the remaining A'^s. Evidently, 
it can be associated with the spectral parameter (e.g. see previous section). 
The Hodge star operator in Eq.(4.6) is needed to make this equation as an 
equation for one-forms The obtained results fit nicely into Cauchy formulation 
of dynamics of both Y-M and gravity. Indeed, under conditions analogous to 
that discussed in [45,46] the space-time (4- manifold) is decomposable into direct 
product M X i? (a trivial fiber bundle) in such a way that all differential op- 
erations acting on 4-manifold are been projected down to 3-manifold M. This 
is essential part of Floer's theory. Furthermore, since SCS{A)/6A — F(A) the 
above Eq.(4.6) can be equivalently rewritten as 

■^A = *[SCS{A)/SA] (4.7) 

so that the Chern-Simons functional is playing a role of a " Hamiltonian" in 
Eq.s(4.7). From the theory of dynamical systems it follows then that the dy- 
namics is taking place between the points of equilibria defined by zero curvature 
condition F(A) = 0. At the same time, using our work, Ref.[50], it is easily rec- 
ognize Eq.(4.7) as an equation for the gradient flow, e.g. see Eq.s(3.14). For the 
sake of space we shall not discuss this topic any further. Interested readers are 
encouraged to consult Ref.[65]. For supersymmetric Y-M fields participating 
in Seiberg-Witten theory the gradient flow equations are discussed in detail in 
Ref.[66] 

The mechanical system described by Eq.(4.7) should be eventually quan- 
tized. Since the quantization procedure is outlined in Ref.[67], to avoid du- 
plications, we shall concentrate attention of our readers on aspects of Floer's 
theory not covered in [67] but still relevant to this paper. To do so, we follow 
Donaldson [11]. This is accomplished in several steps. 

First, in the previous section we noticed that the axially symmetric self-dual 
solution for Y-M fields does not depend on xq (or r) variable. Therefore, if such 
solution is substituted back into Y-M functional, it produces divergent result. 
Although the cure for this issue is discussed in subsection 4.4, in this subsection 
we provide needed background. For this purpose, following Ref. [68] we consider 
the Y-M action S[F] for the pure Y-M ficltO 

S[F] = -\ I d^xtr{F^,F^,). (4.8) 

The duality conditiorF^f — ^e^^apFap allows us then to rewrite this action 
as follows 

^[F] = -4 / d^^HiFf.. T *F^.){.F^. T *F^.)) ± 2tr(F^, * F^,)] (4.9) 

lO JR'I 

Strictly following notations of Ref. [68] we do not indicate that in general the integration 
should be made over some 4-manifold M. In physics literature, and in Eq.s(2.11), it is assumed 
that we are dealing with R'* (or 5'* upon compactification). In Floer's theory it is essential 
that the 4-manifold is decomposable as M y. R . This decomposition should be treated with 
care as described in the Donaldson's book [11] 
^^With the convention that £1234 = —1- 
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since tr{Fi^uF^^) = tr{*F^^u * FtJ.u)- The winding number N for SU(2) gauge 
field is defined aJ^ 

N^-^j d^xtr{F^,*Ff,,) = --^ [ tr{F^,AF^,) (4.10) 

so that use of this definition in Eq.s(4.8),(4.9) produces 

S[F] > TT^ \N\ (4.11) 

with the equahty taking place when the (anti) self-duality condition (e.g. see 
Eq.(2.10) ) holds. In such a case the saddle point action is becoming just a 
winding number (up to a constant). 

Second, if our space-time 4-manifold A4 can be decomposed as M x [0, 1], 
the following identity can be used [11] 

/ tr{Ff,^AF^^)^ I tr(A AdA+^AA AA A) 4= CS'(A). (4.12) 

Jmx[q.i] Jm 3 

Here the symbol = means "up to a constant". The decomposition M x [0, 1] 
reflects the fact that the C-S functional is defined up to mod Z. This ambiguity 
can be removed if we agree to consider C-S functional as a quotient R/Z. Ac- 
cordingly, this allows us to replace M x Rhy M x [0, 1]. Details can be found 
in Ref.[ll]. Thus, one way or another the winding number TV in Eq.(4.10) can 
be replaced by the Chern-Simons functional. 

Third, since the equation of motion for the C-S functional is zero curvature 
condition F = 0, i.e. 

dA + AAA = 0, (4.13) 

implying that the connection A is flat, we can use this result in Eq.(4.12) in 
order to rewrite it as (e.g. for SU(2)) 

If 2 If 

— r / tr(AAdA-|--AA AA A) = --— r / ir(AAAAA). (4.14) 
Stt^ Jj^j 3 247r^ Jj^j 

For other groups the prefactor and the domain of integration will be different 
in general. 

Fourth, zero curvature Eq.(4.13) involves connections which are functions of 
three arguments and a spectral/time parameter. In such setting minimization 
of Y-M functional is not divergent in view of Eq.(4.11). 

Fifth, the obtained result, Eq.(4.14), coincides with that known for the 
winding number for SU(2) instantons in physics literature [8, 19] where it was 
obtained with help of entirely different arguments. It should be noted though 
that in spite of apparent simplicity of these results, actual calculations of C-S 
functionals (invariants) for different 3-manifolds are, in fact, very sophisticated 
[69,70]. In accord with Floer and Ref. [67], we conclude that nonperturbatively 



^We follows notations of Ref. [68] in which is actually standing for £ SU{2) 
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the 4-dimensional pure Y-M quantum field theory is a topological field theory 
of C-S type. 

Sixth, the isomorphism noticed by Louis Witten acquires now natural expla- 
nation. It becomes possible in view of results just presented, on one hand, and 
the fact that AfSn-Elg] = CS{A) (previous section), on another. For fields 
with axial symmetry, equations of motion, Eq.(4.13), for gravity and Y-M fields 
coincide. 

Seventh, the instantons in Floer's theory are not the same as considered 
in physics literature [8,19]. To understand this, we must take into account 
that in Floer's theory manifolds under consideration are 4-manifolds Ai with 
tubular ends. Such manifolds are complete Riemannian manifolds with finite 
number of tubular ends made of half tubes (0, oo) so that locally each such 
manifold looks like Ui = LiX (0, oo) with Li being a compact 3-manifold (called 
a "crossection of a tube" ) and i numbering the tubes. The closure of [J"^i Ui 
is a compact manifold with boundary. If the crossection is S^, then U is 
conformally equivalent to a punctured ball B'^\{0}. This implies that a manifold 
A4 with tubular ends is conformally equivalent to a punctured manifold A4 \ 
{ pi,...,pn} where M is compact. The instanton moduli problem for Ai is 
equivalent to that for the punctured manifold [11]. Recall that the moduli space 
of instantons is defined as set of solutions of anti self-dual equations modulo 
gauge equivalence. 

Being armed with these definitions and taking into account that the (anti) 
self-duality Eq.(4.7) we can interpret the instanton as a path connecting one 
flat connection F = at "time" r = — oo with another flat connection at "time" 
T — oo [11]. It is permissible for the path to begin at one flat connection, to 
wind around a tube (modulo gauge equivalence) and to end up at the same flat 
connection, Ref.[ll], page 22. Evidently, this case involves 4-manifolds with just 
one tubular end. Physically, each flat connection F = represents the vacuum 
state so that the instantons discussed in the Introduction should be connecting 
different vacua. In this sense there is a difference between the interpretation of 
instantons in mathematics and physics literature. As in the case of standard 
quantum mechanics, only imposition of some additional physical constraints 
permits us to select between all possible solutions only those which are physically 
relevant. In the present context it is known that all exactly integrable systems 
are described by the zero curvature equation F = [5,6]. It is also known 
that differences between these equations are caused in part by differences in a 
way the spectral parameter enters into these equations. Since for the Floer's 
instantons F 7^ 0, it means that the curvature F should be parametrized in such 
a way that the "time" parameter should become a spectral parameter when 
F = 0. In this work we do not investigate this problerrF^. Instead, we shall 
focus our attention on different vacua, that is on different (knot-like) solutions 
of zero curvature equation F = Op^ 

^''^See Ref. [7] for introduction into this topic. 

complement of each knot in is 3-manifold. Floer's instantons are in fact connecting 
various three-manifolds. These 3- manifolds (with tubular ends) should be glued together to 
form A^.The gluing procedure is extremely delicate mathematical operation [11]. It is above 
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4.2 The Faddeev-Skyrme model and vacuum states of the 
Y-M functional 



In the light of results just presented, we would like to argue that the F-S model 
is indeed capable of representing the vacuum states of pure Y-M fields. For 
this purpose it is sufficient to recall the key results of the paper by Auckly and 
Kapitansky [71]. These authors were able to rewrite the Faddeev functional 

E[n] = j dv{\dn\^ + \du A dnf} (4.15) 

in the equivalent form given by 

E^[a] = J dv{\D^(pf + \D^(p A D^cpf}. (4.16) 

S3 

In this expression, the covariant derivative -DaV = dip+ [a, ip]. Evidently, E^[a] 
acquires its minimum when (p = a and the connection becomes flat (that is 
covariant derivative becomes zero). Since this result is compatible with those 
discussed in previous subsection, it implies that, indeed, Faddeev's model can 
be used for description of vacuum states for pure Y-M fields. The only question 
remains: Is this model the only model describing QCD vacuum? In view of Eq.s 
(3. 18), (3. 19) it should be clear that this is not the ease. The full explanation 
is given below, in Sections 5,6. In addition, the disadvantage of the F-S 
model as such (that is without modifications) lies in the absence of gap upon 
its quantization as was recognized already by Faddeev and Niemi in Ref.[25]. 
In Sections 5,6 we shall eliminate this deficiency in a way different from that 
described in the Introduction (e.g. in Ref.s[25,26]). In the meantime, we would 
like to find the place for monopoles in our calculations. 



4.3 Monopoles and the Ernst equation 
4.3.1 Monopoles versus instantons 

To introduce notations and for the sake of uninterrupted reading, we need to 
describe briefiy the alternative point of view at the results of previous subsection. 
For this purpose, following Manton [49], we need to make a comparison between 
the Lagrangians for SU(2) Y-M and the Y-M-Higgs fields described respectively 

by 

jCy-m = —triF^^F^^'') (4.17) 

and 

£y-m-h = -\tr{F^,F>'n - ^tr {D^^-D^^^} - ^(1 - * • ^ f (4.18) 

the level of rigor of this paper. To imagine the connected sum of knots [20] is much easier 
than the connected sum of S-manifolds. This sum has physical meaning discussed in Section 
6. 
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with covariant derivative for the Higgs field defined as Z)^$ = + [A^, $] 
and with connection used to define the Y-M curvature tensor F^i, = S^A^ — 
d^K^ + [A^, A^], provided that * =$'^r, A,, = Alt", and = -2eabct''. 

Now the self-duality condition F = *F can be equivalently rewritten as F^j = 
—SijkFko with indices fc running over 1,2,3. Incidentally, in the temporal 
gauge this result is equivalent to Floer's Eq.(4.6) Consider now the limit A — >■ in 
Eq.(4.18). In the Minkowski spacetime the field equations originating from the 
Y-M-Higgs Lagrangian can be solved by using the Bogomolny ansatz equations 
Fij — —EijkDk^ in which Aq = (temporal gauge). Instead of imposing the 
temporal gauge condition, we can identify the Higgs field $ with Ao so that 
the Bogomolny equations read now as follows: 

F„ = -Eyfci^fcAo. (4.19) 

Bogomolny demonstrated that the Prasad-Sommerfield monopole solution can 
be obtained using Eq.(4.19). Thus, any static (that is time- independent) so- 
lution of self-duality equations is leading to Bogomolny-Prasad-Sommerfield 
(BPS) monopole solution of the Y-M fields, provided that we interpret the 
component Aq as the Higgs field. Suppose now that there is an axial symme- 
try. Forgacs, Horvath and Palla [72] (FHP) demonstrated equivalence of the 
set of axially symmetric Bogomolny Eq.s (4.19) to the Ernst equation. The 
static monopole solution is time-independent self-dual gauge field. Because of 
this time independence, its four-dimensional action is infinite (because of the 
time translational invariance) while that for instantons is Enite. Furthermore, 
the boundary conditions for monopoles and instantons are different. The infin- 
ity problem for monopoles can be cured somehow by considering the monopole 
dynamics [68] but this topic at this moment "is more art than science", e.g. 
read [68], page 309. For the same reason we avoid in this section talking about 
dyons (pseudo particles having both electric and magnetic charge). Hence, we 
would like to conclude our discussion with description of more mathematically 
rigorous treatments. By doing so we shall establish connections with results 
presented in previous sections 

4.4 Calorons. 

Calorons are instantons on x S^. From this definition it follows that, phys- 
ically, these are just instantons at finite temperaturj^. Calorons are related 
to both instantons on R'' (or S*^) and monopoles on R'^(or S^). Heuristically, 
the large period calorons are instantons while the small period calorons are 
monopoles [73,74]. These results do not account yet for the fact that both the 
Y-M action and the self-duality equations are conformally invariant. Atiyah 
[75] . noticed that the Euclidean metric can be represented either as 

dsl = {dx^f + [dx^'f + [dx^f + [dx^f (4.20a) 

^^This explains the word "caloron". 
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or as 

ds^ = -^[R^i- + R^dif^] (4.20b) 

with R being some constant. The above representation involves polar r,(p 
coordinates in the {x^,x'^) plane thus implying some kind of axial symmetry. 

2 

Since self-duality equations are conformally invariant, the prefactor can be 
dropped so that the Euclidean space R'' becomes conformally equivalent to 
the product H"^ x S'^.For such manifold the constant scalar curvature of the 
hyperbolic 3-space H"^ is —1/B?. Furthermore, the remaining term represents 
the metric on a circle of radius R. Following Ref.[73], let (x^, x^, x'^) be coor- 
dinates for the hyperbolic ball model of H'^ so that the radial coordinate be 

7^ = \J{x'^f + {x'^f + {x'^f. Let < 7^ < i?. Let r be a coordinate on 
with period then the metric on H"^ x can be represented as 

ds\j = dr^ + A2(rf7^2 + n^dn^), (4.21a) 

where A = (1 — Tl/R)~^ and dfi^ is the metric on 2-dimensional sphere. If 
we introduce an auxiliary coordinate /i = (i?/2) arctanh(7?./i?), and complex 
coordinate z = iJ, + iT,the above hyperbolic metric can be rewritten as 

ds]j = dr^ + djji^ + S^dfi^ (4.21b) 

with S = (ii/2) sinh(2/x/i?). By analogy with transition from Eq.(4.20a) to 

(4.20b) we can proceed as follows. Let r = ^ {y^)^ + (y^)^ + (y^)^ with {y^, y'^, y^, 
being coordinates on R^. By letting t = y° the Euclidean metric can be written 
as usual, i.e. 

ds% = dt^ + dr^ + r'^d^f (4.22a) 

so that 

dsjj = ^^dsl, (4.22b) 

with ^ = (E/2)[cosh(2/x/i?) + cos(2T/i?)]. This correspondence between R^'xR^ 
and x is made with help of the mapping w = tanh{z/R) (with w = r + it 
and 13 = ttR) . Let M xS^ (or x R) then, in view of conformal invariance, 
we can rewrite Eq.(4.8) as 

S[F] = -l [ tr{F^^F^''')E'^dTdiJ,dn. (4.23) 

We have to rewrite the winding number, Eq.(4.10), accordingly. Since it is 
a topological invariant, this means that the self-duality equations must be ad- 
justed accordingly. For instance, for the hyperbolic calorons on x the 
self-duality equation reads 

Foi = ^SijkFjk- (4.24) 

The action S'[F] now is finite with tr{Fni,F'^'') — >■ when — >■ oo. For hyperbolic 
instantons we have finite action with tr{FfivF^^'^) — >■ when /x^ + — >■ oo. 
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The results just described match nicely with the results by Witten [76] on 
Euclidean SU(2) instantons invariant under the action of S0(3) Lie group. His 
results will be discussed in detail in the next section. Notice, that Euclidean 
metric, Eq.(4.22a), becomes that for x if we rewrite it as 

dsl = r\^^^^ + (4.25a) 

and, as before, we drop the conformal factor so that it becomes 

rf4 = ^*'+/''' +dn^. (4.25b) 



Interestingly enough, that results by Witten initially developed for ¥P x 5^ 
can be also used without change for x S*^ and H"^ x R since the action of 
S0(3) pulls back to these manifolds [73] . This fact is of importance since such an 
extension makes his results compatible with both Floor's method of calculations 
for Y-M fields and with results of Section 3. Omitting all details, the action, 
Eq.(4.23), is reduced to that known for two dimensional Abelian Ginzburg- 
Landau (G-L) model "living" on the hyperbolic 2 manifold X coordinatized by 
H and T with the metric 

dsl = . (4.26) 

Explicitly, such G-L action functional Sg-l is given by [73] 

Sg-l = \j dTdi^[E'{V X A)2 + 2 |(V + iA)4>\' + - |0|')] (4.27) 

X 

with A and (p being respectively the Abelian gauge and the Higgs fields, (p = 
(j)i+i<p2 so that \4>\^ = (j)1 + (j)2- This functional is obtained upon substitution of 
solution of the self-duality equations into the Y-M action functional, Eq.(4.23). 
We refer our readers to the original paper, Ref.[73] for details. In the limit 
/3 — >■ oo the above hmctional coincides with that obtained by Witten [76]. The 
self-duality equations obtained by Witten describe instantons which lie along the 
fixed axis while Fairlie, Corrigan, 't Hooft, and Wilczek [77] developed an ansatz 
(CFtHW ansatz) for the self-duality equations producing instantons at arbitrary 
locations. Manton [78] demonstrated that Witten's and GFtHW multi-instanton 
solutions are gauge equivalent while Harland [73,74] demonstrated how these 
instantons and monopoles can be obtained from calorons in various limits. The 
obtained results provide needed background information for solution of the gap 
problem. This solution is discussed in the next section. 



5 Solution of the gap problem 

5.1 Idea of the proof 

By cleverly using symmetry of the problem Witten [76] reduced the non Abelian 
Y-M action functional to that for the Abelian G-L model "living" in the hyper- 
bolic plane. This is one of examples of the Abelian reduction of QCD discussed 
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in our paper, Ref.[79]. Vortices existing in the G-L model could be visualized 
as made of some two-dimensional surfaces (closed strings) living in the ambi- 
ent space-time. These are known as Nambu-Gotto strings. Their treatment 
by Polyakov [80] made them to exist in spaces of higher dimensionality. In or- 
der for them to be useful for QCD, Polyakov suggested to modify string action 
by adding an extra (rigidity) term into string action functional. By doing so 
the problem was created of reproducing Polyakov rigid string model from QCD 
action functional. The latest proposal by Polyakov [81] involves consideration 
of spin chain models while that by Kondo [82] involves the F-S model derived 
directly from QCD action functional. As explained in[79], in the case of scatter- 
ing processes of high energy physics one is confronted essentially with the same 
combinatorial problems as were encountered at the birth of quantum mechanics. 
In Ref.[83] we explained in detail why Heisenberg's (combinatorial) method of 
developing quantum mechanical formalism is superior to that by Schrodinger. 
In Ref.[74] using these general results we demonstrated how the combinatorial 
analysis of scattering data leads to spin chain models as microscopic models 
describing excitation spectrum of QCD. Thus, the mass gap problem can be 
considered as already solved in principle. Nevertheless, in [94] such a solution 
is obtained "externally", just based on the rules of combinatorics. As with 
quantum mechanics, where atomic model is used to test Heisenberg's ideas, 
there is a need to reproduce this combinatorial result " internally" by using mi- 
croscopic model of QCD. For this purpose, we shall use the G-L functional, 
Eq.(4.27). By analogy with the flat case, we expect that it can be rewritten in 
terms of interacting vortices. In the present case, in view of Eq.(4.26), vortices 
"live" not in the Euclidean plane but in 3-1-1 Minkowski space-time. This is 
easy to understand if we recall the SO(3)<=^SU(2) correspondence and take into 
account the analogous correspondence between SU(1,1) and S0(2,l). 

Within such a picture it is sufficient to look at evolution dynamics of the 
individual vortex. Typically, it is well described by the dynamics of the continu- 
ous Heisenberg spin chain model [84,85] in Euclidean space. In the present case, 
this formalism should be extended to the Minkowski space and, eventually, to 
hyperbolic space (that is to the case of Abelian model discovered by Witten). 
Details of such an extension are summarized in Appendix B. After that, the 
next task lies in connecting these results with the Ernst equation. In the next 
subsection we initiate this study. 

5.2 Heisenberg spin chain model and the Ernst equation 

For the sake of space, this subsection is written under assumption that our read- 
ers are familiar with the book " Hamiltonian methods in the theory of solitons" 
[86] (or its equivalent) where all needed details can be found. The contin- 
uous XXX Heisenberg spin chain is described with help of the spin vectoiF^I 

^"In compliance with [86] we suppress the time-dependence. 
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S{x) = {Si{x), S2{x), Ss{x)) restricted to live on the unit sphere 5^ 



S\x)=J2s!{x) = l (5.1) 

i=l 

while obeying the equation of motion 

d ^ - ^2 ^ 

known as the Landau-Lifshitz (L-L) equation. By introducing matrices U{X) 
and V{X) via 

U{X) = ^5, V{X) = ^5 + ^S-^S, S = S-a (5.3) 

so that ai is one of Pauli's spin matrices and A is the spectral parameter and 
requiring that = I, where / is the unit matrix, the zero curvature condition 

l-U-lv+[U,V]=, (5.4) 

is obtained. With account of the constraint 5^ = / it can be converted into 
equation 

equivalent to Eq.(5.2). The correspondence between Eq.s(5.2) and (5.5) can be 

made for S{x, t) matrices of arbitrary dimension. 

Having in mind Witten's result [76] , we want now to extend these Euclidean 
results to the case of noncompact Heisenberg spin chain model "living" either 
in Minkowski or hyperbolic space. In doing so wc follow, in part, Rcf.[56] and 
Appendix B. For this purpose we need to remind our readers some facts about 
the Lie group SU(1,1). Since this group is related to S0(2,l), very much like 
SU(2) is related to S0(3), we can proceed by employing the noticed analogy. 
In particular, since S = S ■ a, we can preserve this relation by writing now 
S = S ■ T. Using this result we obtain, 

^ = ( fs+ -iS^ ) ^ = ± (5-6) 

where the form of matrices generating su(l,l) Lie algebra is similar to that 
for Pauli matrices. This time, however, dot S* = —1 even though = I . 
Explicitly, {S^)^ — (5^)^ — (5^)^ = l,that is the motion is taking place on the 
unit pseudosphere 6"^'^. Matrices Tj generating su(l, 1) are fully characterized 
by the following two properties 

tr{TaTp) =2gai3 , [tcT/j] = 2ifai3jTj; ga0 = diag{-l,-l,l);a,P,j = 1,2,3 

(5.7) 
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with fap^ being structure constants for su{l, 1) algebra. An analog of the 
equation of motion, Eq.(5.5), now reads 

po = J2r''S,^S,. (5.8) 

/3,7 

If we defines the Poisson brackets as {S'"(x), 5'^(y)} = - f^^f S'l {x)6{x - y), 
then the above equation of motion can be rewritten in the Hamiltonian form 

^S^ = {H,S"}, (5.9) 
provided that the Hamiltonian H is given by 

oo oo 

H=l J dx{V,S")gc0{V,S^) = ^tr J dx{V,Sf . (5.10) 

— oo — oo 

Since now the motion takes place on pseudosphere S"^, it is convenient to intro- 
duce the pseudospherical coordinates by analogy with spherical, e.g. 

S^{x,t) = sinh 6{x,t)cosip{x,t),S^{x,t) = smh6{x,t)sirnp{x,t),S''{x,t) = cosh6{x,t). 

(5.11) 

Also, by analogy with spherical case we can use the stereographic projection 
: from pseudosphere to hyperbolic plane. Recall [102], that in the case of a 
sphere S*^ the inverse stereographic projection: from complex plane C to is 
given by 

The mapping from C to is obtained with help of Eq.(5.12) in a straightfor- 
ward way as 
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— n t S — o^S — o* (5.13) 

i-i^r i-iei' i-i^i' 

Using this correspondence the equations of motion, Eq.(5.10), rewritten in 
terms of ^ and 5* variables (while keeping in mind that they are parametrized 
by X and t) are given by 



'de'^ dx^^~^ 1- \dx' 



In the static (f— independent) case the above equation is reduced to 

del' - l)V^e = 2r (V.O' (5.15) 

easily recognizable as the Ernst equation. In his paper, Ref. [42], Ernst used 
variational principle applied to the functional Eq.(3.18). From Appendix B we 
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know that both the L-L equation and its hypcrbohc version describe the motion 
of (could be knotted) vortex filament. Because of this, the functional, Eq.(3.18), 
should undergo the same reduction as was made in going from Eq.(Bl.a) to 
(Bl.b). Explicitly, this means that the functional, Eq.(3.18), should be reduced 
in such a way that the Hamiltonian, Eq.(5.10), should be replaced by 

H = -2 f dx ^^''^t . (5-16) 
J (1-1^1^)2 

— oo 

where the sign in front is chosen in accord with Ref. [87] and our Eq.(3.22). 
The Hamiltonian equation of motion, Eq.(5.9), reproducing Eq.(5.14) can be 
obtained if the Poisson bracket is defined as by {^{x), Ciu)} = (1 ~ ICI^)^'5(2; — 
y).The obtained results set up the stage for quantization. It will be discussed in 
subsection 5.4. In the meantime, we need to connect results of Witten's work, 
Ref. [76] , with those we just obtained. 



5.3 From Abelian Higgs to Heisenberg spin chain model 
5.3.1 The Abelicin Higgs model 

The work by Witten [76] had been further analyzed in the paper by Forgacs 
and Manton [88]. The major outcome of their work lies in demonstration 
of uniqueness of the self-duality ansatz proposed by Witten. The self-duality 

equations obtained in Witten's work are reduced to the system of three coupled 
equations describing interaction between the Abelian Y-M and Higgs fields 

doipi + Ao(p2 = di'P2 - ^iViJ (5.17a) 

diipi + Aiip2^ -ido(p2- Aoipi), (5.17b) 

r\doAi - diAo) = 1 - - (5.17c) 

To analyze these equations, we recall that the original self-duality equations for 
the Y-M fields are conformally invariant. We can take advantage of this fac;t now 
by temporarily dropping the conformal factor in Eq.(5.17c). Then, the above 
equations become the Bogomolny equations for the flat space Abelian Higgs 
model, e.g. for the model described by the action functional. Eq.(4.24), with the 
conformal factor S = 1 [89]. Such obtained equations contain all information 
about the Abelian Higgs model and, hence, they are equivalent to this model. 
It is of importance for us to demonstrate this explicitly for both Euclidean 
and hyperbolic spaces. For this purpose we introduce a covariant derivative 
Dfi = df_i — iA^, /X = 0, 1, and the complex field 4> = (ji^ + i(j)2- Consider the 
Bogomolny equation following [68]: 

Do(l> + iDicl) = 0. (5.18) 

Using the above definitions straightforward computation reproduces Eq.s(5.17a,b) 
These equations can be used to obtain 

{Do - iDi) {Do + iDi)(^ = (5.19) 
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implying 



(5.20) 

where the last equality was obtained with help of Eq.(5.17c). Evidently, the 
equation 

(DoDo + DiDi)<j) +1(1-^1- ^l)4> ^ (5.21) 

is one of the equations of " motion" for the G-L model on , e.g. see Ref.[89](Eq.(11.3) 
page 98). The second is the Ampere's equation 

et,,d^,{r''B) - i{^D,~^ - 0D,0) (5.22) 

with the "magnetic field" B — d^Ai — diA^. Details of derivation are given in 
Ref.[68], pages 198-199. Eq.(5.22) also coincides with that given in the book by 
Taubs and Jaffe, Ref.[105] (Eq.(11.3) page 98). 

Corollary 1. Since both equations can be obtained by mimization of the 
functional, Eq.(4.27), they are equivalent to the Abelian Higgs model which, in 
turn, is the reduced form of the Y-M functional for pure gauge fields. 

We continue with the discussion of Witten's treatment of Eq.s(5.17) since 
we shall need his results later on. First, he selects physically convenient gauge 
condition via df^A^^ = 0. This leads to the choice: — Sfj^^d^ip (for some scalar 
ijj)- With such a choice for A^^ the first two of Eq.s(5.17) can be rewritten as 

{do - dotp)(pi = {di - ditp)(f2, (5.23a) 

{di - di%l))ip^ = -{do - do^p)(p2- (5.23b) 

Let now = e^Xi ^iid = e'^X2- Then the above equations are reduced to 
the Cauchy-Riemann-type equations: d^Xi = diX2 ^^'^ ^^iXi — (^0X2- Introduce 
the function f = Xi ~ *X2- Then, the last of Eq.s(5.17) acquires the form 

-r^vV^ l-//*e2'^. (5.24) 

Notice that -r^V^ = -r^(^ + is the hyperbohc Laplacian [90]. Eq.(5.24) 
is still gauge invariant in the sense that by changing / — fh and ip ^ tp — 
5 ln{hh* ) in this equation we observe that it preserves its original form. This 
is so because V^ln(/i/i*) = for any analytic function which does not have 
zeros. If h does have zeros for r > 0, then substitution oi ^ ip — ^\n{hh*) 
into Eq.(5.24) produces isolated singularities at these zeros. After these remarks, 
Eq.(5.24) can be simplified further. For this purpose, let ip = Inr— ^ ln(//*)-|-p, 
provided that V'^ln(//*) = for any analytic function / which does not have 
zeroj^. Under such conditions we end up with the Liouville equation 

W^p^e^P. (5.25) 

It is of major importance for what follows. 

^^In the case if it does, the treatment is also possible as explained by Witten. Following 
his work, we shall temporarily ignore this option. 
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5.3.2 The Heisenberg spin chain model 

The results of Appendix B imply that the L-L Eq.(5.2) (or their hyperbolic 
equivalent, Eq.(5.8)) could be interpreted in terms of equations for the Serret- 
Prenet moving triad. Treatment along these lines suitable for immediate appli- 
cations is given in papers by Lec and Pashacv [91] and Pashacv [92]. Below we 
superimpose their results with those of our work, Ref. [84] , to achieve our goals. 

We begin with definitions. A collection of smooth vector fields nju(x,t), 
/X = — 2, forming an orthogonal basis is called the "moving frame". If x G 5 
where S is some two dimensional surface, then let ni(x,t) and n2(x,t) form 
basis for the tangent plane to 5 Va; e S. Then, the Gauss map (that is the 
map from S to two dimensional sphere S'^ or pscudosphcrc S^'^) is given by 
n2(x,t)= s. By design, it should obey Eq.(5.1). This observation provides 
needed link between the spin and the moving frame vectors. Details are given 
in [91,92] and Appendix B. It should be clear that since one can draw curves 
on surfaces both formalisms should involve the same elements. The restriction 
for the curve to lie at the surface causes additional complications in general but 
nonessential in the present case. 

Next, we introduce the combinations n± = ng ±mi possessing the following 
properties 

(n+,n+) = (n_,n_) = , (n+,n_) = 2//^^^ (5.26) 

where = 1 for S*^ and = —1 for S^'^ and H^. Furthermore, (.., ..) defines 
the scalar product (in Euclidean or pseudo-Euclidean spaces). Also, 

n+ X s =zn+,n_ x s = — in_,n_ x n+=2iK^s. (5-27) 

In addition, we shall use the vectors 

= y(5/xs,n+) and = y(a^s,n_) (5.28) 

in terms of which the equations of motion for the moving frame vectors look as 
follows: 

D^n+ = -2K^q^s, (5.29) 

a^s =qi^n_ + g^n+, (5.30) 

with covariant derivative = d^ — and = — 2K^(ni, 9^no). Consider 
now Eq.(5.30) for /x = 1. Apply to it the operator di and use the equations of 
motion and the definitions just introduced in order to obtain 

dls = (Diqi) n_ + (Digi) n+ - 4 l<?i l" s. (5.31) 

It can be shown that go = iD\qi. In view of this, Eq.(5.30) for = acquires 
the following form: 

Oqs =iDiqin- — iDiqin^. (5.32) 

This equation happens to be of major importance because of the following. 
Multiply (from the left) Eq.(5.31) by sx and use Eq.s(5.27). Then (depending 
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on signature of k^) the obtained result is equivalent to the L-L Eq.(5.2) or its 
pseudoeuclidean version, Eq.(5.8). Furthermore, for this to happen the fields 
and must be subject to the following constraint equations obtainable directly 
from Eq.s (5.29) 

D^j^qu = D^q^^, (5.33a) 

[Df^.Di,] = -2K^(q^,qi, - qi,qf_,). (5.33b) 

We are going to demonstrate now that these equations are equivalent to Eq.s(5.17) 
obtained by Witten. 

We begin with the following observation. Let indices /i and v be respectively 
1 and 0. Then, taking into account that go ~ iDiqi we can rewrite Eq.(5.33b) 
as 

i^io = Bi = -2Kh{qiDiqi - qiDiqi). (5.34) 

Surely, by symmetry we could use as well: qi — — iDo'Zo- This would give us 
an equation similar to Eq.(5.34). Take now the case = — 1 (that is consider 
S^'^) in these equations and compare them with the Ampere's law, Eq.(5.22). 
We notice that these equations are not the same. However, since the G-L model 
was originally designed for phenomenological (thermodynamical) description of 
superconductivity (as explained in detail in our work, Ref.[84]), we know that the 
underlying equations (obtainable from the G-L functional) contain the London 
equation which readio 

V X B = CB (5.35) 

with C being some constant (determined by physical considerations) . Evidently, 
in view of the London (5.35), Eq.s(5.22) and (5.34) become equivalent. Gonsider 
now Eq. (5.33a). To understand better this equation, it is useful to rewrite 
Eq.(5.18) as follows 

Dq<I) — —iDi<j) or Dqc^i ~ Di4>q, (5.36) 

where ipi = (j) and 0q = Take into account now that (p ~ a + ib and 

identify (f)-^ with qi and 4>q with q^. Then, Eq.(5.33b) acquires the following 
form (k^ = —1) : 

{doVi - diVo) = -I'iiMi - Ml) = -4(a' + (5.37) 

Looking at Eq.(5.17c) we can make the following identifications: Vi — Ai,Vo = 
Ao,±2a — (pi,±2b = Lp2- Then, comparison between Eq.s(5.17c) and (5.37) 
indicates that we are still missing a factor of in the l.h.s. and 1 in the r.h.s. 
Looking at Witten's derivation of the Liouville Eq.(5.25), we realize that these 
two factors are interdependent. By clever choice of the function ip they can 
be made to disappear. This makes physical sense since locally the underlying 
surface is almost flat. This observation makes Eq.s(5.37) and (5.24) (or 5.17c) 
equivalent. 

^^This is not the form of the London equation one can find in textbooks. But in our work, 
Ref. [84], we demonstrated that Eq. (5.35) is equivalent to the London equation. 
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Corollary 2. The L-L and 2 dimensional G-L models are essentially equiv- 
alent in the sense just described both in Euclidean and in Minkowski spaces. 

Corollary 3. The "hyperbole" L-L Eq. {5.1 A) or its Euclidean analog should 
be identified with Floer's Eq.(4.6). 

These results play an important role in the rest of this work and, in partic- 
ular, in the next subsection. 

5.4 The proof (implementation) 
5.4.1 General remarks 

In Ref.[79], we demonstrated how treatment of combinatorial data associated 
with real scattering experiments leads to restoration of the underlying quantum 
mechanical model reproducing the meson spectrum. It was established that 
the underlying microscopic model is the Richardson-Gaudin (R-G) XXX spin 
chain model originally designed for description of spectrum of excitations in the 
Bardeen-Cooper-Schriefer (BCS) model of superconductivity. Subsequently, 
the same model was used for description of spectra of the atomic nuclei. Since 
the energy spectrum of the BCS model has the famous gap between the ground 
and the first excited state, the problem emerges : 

Can spectral properties of nonperturbative quantum Y-M field 
theory be described by the r-g model? 

To answer this question afRrmatively the "equivalence principle" discovered 
by L.Witten is very helpful. Using it, we can proceed with quantization of pure 
Y-M fields by using results by Korotkin and Nicolai, Ref. [31], for gravity. By 
comparing the main results of our paper, Ref. [79], done for QCD, with those of 
Ref. [31], done for gravity, we found a complete agreement. In particular, the 
Knizhnik-Zamolodchikov Eq.s(4.14),(4.15) and the R-G Eq.(4.29) of Ref.[79] 
coincide respectively with Eq.s(4.27),(4.26) and (4.50) of Ref.[31] even though 
methods of deriving of these equations are entirely different! Both Ref.s [79] 
and [31] do not reveal the underlying physics sufficiently deeply though. In the 
remainder of this section we shall explain why this is indeed so and demonstrate 
ways this deficiency can be corrected. Experimentally the challenge lies in de- 
signing scattering experiments providing clean information about the spectrum 
of glueballs. Thus far this task was accomplished only in lattice calculations 
done for unphysically large number of colors, e.g. Nc — cxi.[23]. When it comes 
to interpreting real experiments (always having only three colors to consideiF^) . 
the situation is even less clear, e.g. see Ref. [93]. Hence, the gap problem is full 
of challenges for both theory and experiment. Fortunately, at least theoretically, 
the problem does admit physically meaningful solution as we explained already. 
We continue with ramifications in the next subsection. 

^''E.g. read Section 6. 
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5.4.2 From Landau-Lifshitz to Gross-Pitaevskii equation via Hashimoto 

map 



Since the F-S model is believed to be capable of describing QCD vacua and 
is also capable of describing knotted/linked structures [17], two questions arise: 
a) Is this the only model capable of describing QCD vacua? b) To what extent 
it matters that the F-S model is also capable of describing knots and links? 
The negative answer to the first question follows from Corollary 3 implying 
that, in principle, both Euclidean and hyperbolic versions of the L-L equation 
are capable of describing QCD vacua: different vacua correspond to different 
steady-state solutions of the L-L equations. The negative answer to the second 
question can be found in a review, Rcf.[85], by Annalisa Calini. From this 
reference it follows that, besides the F-S model, knotted/linked structures can 
be also obtained by using standard (that is Euclidean) L-L equation, e.g. see 
Eq.(B.4) of Appendix B. This fact still docs not explain why knots/links are of 
importance to QCD. We address the above issues in more detail in Section 6. In 
view of what is said above, wether or not the hyperbolic version of L-L equation 
is capable of describing knotted structures is not immediately important for 
us. Far more important is the connection between the hyperbolic L-L and the 
Ernst equation. Only with this connection it is possible to reproduce results by 
Korotkin and Nicolai [31]. 

Eq.(3.19) is just the F-S functional without winding number term. When the 
commutation relations for su(l,l) introduced in subsection 5.2 are replaced by 
those for su(2) this leads to the standard L-L equation (instead of Eq.(5.14)). 
This replacement causes us to abandon the connection with Ernst equation 
and, ultimately, with the results of Ref.[31]. In such a case the gap problem 
should be investigated from scratch. In Ref.[25] Faddccv and Niemi indicated 
that, unless some amendments to the F-S model are made, it is gapless. At the 
same time from Appendix B it is known that the L-L equation associated with 
the F-S model can be transformed into the NLSE with help of the Hashimoto 
map. Recently, Ding [94] and Ding and Inoguchi [95] were able to find analogs 
of the Hashimoto map for the vortex filaments in hyperbolic, de Sitter and 
anti de Sitter spaces. It is helpful to describe their findings using terminology 
familiar from physics literature [96]. This leads us to the discussion of properties 
of the Gross-Pitaevskii equation known in mathematics as the NLSE. In the 
system of units in which H = 1 and m = 1/2 this equation can be written as 



Zakharov and Shabat [97,98] performed detailed investigation of this equation 
for both positive and negative values of the coupling constant k. For «; < 
the above equation is used for description of knots/links [85]. The standard 
Hashimoto map brings the L-L equation associated with the truncated F-S 
model to the NLSE with k < [94, 95]. Prom the same references it can be 
found that the Hashimoto-like map brings the (hyperbolic) L-L-like equation to 
the NLSE for which k > 0. Zakharov and Shabat studied in detail differences 



[86] 




(5.38) 
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in treatments of the NLSE for both negative and positive couphng constants. 
This difference is caused by difference in underlying physics which in both cases 
can be explained in terms of the properties of non ideal Bose gas [99,100]. The 
attentive reader might have noticed at this point that Eq.(5.38) apparently 
contains no information about the number of particles in such a gas. This 
parameter, in fact, is hidden in the constant c (the chemical potential) or it can 
be obtained selfconsistently with help of Eq.(5.38) (from which c is removed in a 
way described in Appendix B) as explained in Ref. [100]. With this information 
at our disposal we are ready to make the next step. 

5.4.3 From non ideal Bose gas to Richardson-Gaudin equations 

Even though statistical mechanics of 1-d interacting Bose gas was considered in 
detail by Lieb and Linger [101], solid state physics literature is full of refinements 
of their results up to moment of this writing. These refinements have been 
inspired by experimental and theoretical advancements in the theory of Bose 
condensation [96]. Among this literature we selected Ref.s[102,103] as the most 
relevant to our needs. 

Following [102], the Hamiltonian for N interacting bosons moving on the 
circle of length L is given by 



with constant 2c coinciding with 2k in the system of units H = 1 and m = 1/2. 
The case c > (repulsive Bose gas) corresponding to the L-L equation in the 
hyperbolic plane/space happens to be of immediate relevance. Only for this case 
it is possible to establish the connection with work by Korotkin and Nicolai [31]! 

We begin by noticing that in the standard BCS theory of superconductivity 
electrons are paired into singlets (Cooper pairs) with zero centre of mass mo- 
mentum. The pairing interaction term in this theory accounts only for pairs 
of attractive electrons with opposite spin and momenta so that the degener- 
acy for each energy state is a doublet, with level degeneracy D, — 2 In the 
interacting repulsive Bose gas model byRichardson [104] to mimic this pairing 
he coupled two bosons with opposite momenta ±kj into one (pseudo) Cooper 
pair. An assembly of such formed pairs forms repulsive Bose gas which in the 
simplest case is described by the Hamiltonian, Eq.(5.39). Hence, the fermionic 
BCS-type model with strong attractive pairing interaction can be mapped into 
bosonic repulsive model proposed by Richardson. Although the idea of such 
mapping looks very convincing, its actual implementation in Ref. [102] has some 
flaws. Because of this, we shall use results of this reference selectively. For this 
purpose, fist of all we need to make an explicit connection between the repulsive 
Bose gas model described by Eq.(5.39) and the model proposed by Richardson. 
In the weak coupling limit cL 1 the Bethe ansatz equations for the repulsive 

■^^We use here the same notations as in our work, Ref. [94]. 




(5.39) 
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Bose gas model described by the Hamiltonian, Eq.(5.39), acquire the following 
form: 

N 

fc. = -^ + -^^,.= l,...,Ar. (5.40) 

Here = 0, ±1, ±2, ... denote the excited states for fixed N. To link this result 
with Richardson's (repulsive boson) model, consider the case of even number of 
bosons and make N = 2M. Next, consider the ground state of this model first. 
To the first order in c, it is clear that we can write ki = ±y/Ei. Specifically, 
let fci,2 = ±^/Ei, ^3,4 = ±\/E3, k2M-i,2M = ±\/Em . Using these results in 
Eq.(5.40), with the accuracy just stated, the Bethe ansatz equations after some 
algebra are converted into the following form: 



2c ^ - Ei 2Ei 

j — 1 



, M; M < M. (5.41) 



To analyze these equations, we expect that our readers are familiar with works 
of both Richardson-Sherman, Ref.[105], and Richardson, Ref.[104]. In [105] 
diagonalization of the pairing force Hamiltonian describing the BCS-type su- 
perconductivity was made. Such a Hamiltonian is given by 

H = Y, - 5^' b\bf,, (5.42) 

/ 

where Nf = i(aj^aj_+a^_a/_), bf ~ aj^aj^, with a|-^ and ajcr being fcrmion 
creation and annihilation operators obeying usual anticommutation relations 
[a/CT) = S^„'dff', where tr = ± denotes states conjugate under time 
reversal. The above Hamiltonian is diagonalizcd along with the seniority oper- 
ators (taking care of the number of unpaired fermions at each level /) defined 

by 

!>/ = 4_^a/- - 4_a/- . (5.43) 

By construction, [H,Nf] = [H,i>f] = 0. The classification of the energy levels 
is done in such a way that the eigenvalues Vf of the operator 1^/(0 and a) are 
appropriate for the case when g = 0. This observation allows us to subdivide the 
Hamiltonian into two parts: i.e. that which does not contain Cooper pairs 
(for which i/f = a) and i?2,«-e.that which may contain such pairs (for which 
vf =0). The matrix elements of H2 are calculated with help of the bosonic-type 
commutation relations 

[br,Nf'] = Sff'bf and [bf , b\,] = Sfril - 2Nf,). (5.44) 

These commutators are bosonic but nontraditional. In the traditional case we 
have [6/ , = Sjf. We refer our readers to Ref.[105] for details of how this 
commutator difficulty is resolved. In the light of this resolution, in Ref.[104] 
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Richardson proposed to deal with the interacting bosons model from the be- 
ginning. Supposedly, such bosonic model can be designed to reproduce results 
of the fermionic pairing model of Ref. [105]. An attempt to do just this was 
made in Ref. [102]. In the repulsive boson model by Richardson the "pairing" 
Hamiltonian is given b'vF^ 

H = Y^ 2em + f E^, 4^/'- (5.45) 

in which hi and Afi are bosonic analogs of Nf and bf. It is essential that 
the sign of the coupling constant g is nonnegative (repulsive bosons). Upon 
diagonalization, the total energy E is given by 

n m 

E^Y.^m + Y.^, (5.46) 

1=1 j=i 

so that summation in the first sum takes place over the unpaired bosons while 
in the second- over the paired bosons whose energies Ej are determined from 
the Richardson's equation (Eq.(2.29) of Ref. [104] 
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2g ^ 2ei ~ Ek E, ~ Ek 



E 

i=l 



0,fc = l,...,TO (5.47) 



in which n is the total number of single particle (unpaired) levels, m is the total 
number of pairs, di = ^{2^1 + ^i)- From [104] it follows that for the bosonic 
model to mimic the BCS-type pairing model the degeneracy factor ili — 1 and 
1^1 — 0. It should be noted though that such an identification is not of much 
help in comparing the repulsive bosonic model with the attractive BCS-type 
fermionic model (contrary to claims made in Ref. [102]). This can be easily 
seen by comparison between Eq.(5.47) (that is Eq.(2.29) of Ref. [104]) with such 
chosen fli and with Eq.(3.24) of Ref. [105]. By replacing g in Eq.(5.47) by 
—g we still will not obtain the analog of the key Eq.(3.24) of Ref.[105]! This 
fact has group-theoretic origin to be explained in the next subsection. In the 
meantime, Eq.(5.47) still can be used to connect it with Eq.(5.41) originating 
from different bosonic model described by the Hamiltonian Eq.(5.39). To do so 
we follow the path different from that suggested in Ref. [102]. Instead, following 
the original Richardson's paper [104], we let n = 1 in Eq.(5.47) then, without 
loosing generality, we can put ei = so that Eq.(5.47) acquires the following 
form 

1 1 2 

2 — 1 



^^To avoid ambiguities, our coupling constant ^ is chosen exactly the same as in [104], 
^^Since Gaudin's equation is obtained in the limit \g\ oo. from Eq.(5.47) the spin -like 
model described by this equation is known as the Richardson-Gaudin (R-G) model. 
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The rationale for replacing to by M is given on page 1334 of [104]. Evidently, 
Eq.s (5.41) and (5.48) are identical. This observation allows us to use the 
Richardson model instead of that described by Eq.(5.39). At first sight such 
an identification looks a bit artificial. To convince our readers that it does 
make sense, we would like to use the work by Dhar and Shastry [106,107] on 
excitation spectrum of the ferromagnetic Heisenberg spin chain. By analogy 
with Eq.(5.41) these authors derived a similar equation obtained by reducing 
the Bethe ansatz equations for Heisenberg ferromagnetic chain. It readj^ 

2 — 1 

Even though Eq.s(5.48) and (5.49) look almost the same, they are not the same! 
The crucial difference lies in the signs in front of the second term in the r.h.s. of 
these equations. Because of this difference Heisenberg's ferromagnetic spin chain 
model is mapped onto Bose gas model with attractive interaction in complete 
accord with what was said immediately after Eq.(5.38). Regrettably, this result 
is still not the same as for the BCS-type model investigated in Richardson- 
Sherman's paper, Ref.[105]. This fact was recognized and discussed in some 
detail already by Richardson [104]. For completeness, we mention that the 
problem of BCS-Bose-Einstein condensation (BEC) crossover which follows 
exactly the qualitative picture just described was made quantitative only 
very recently in Ref.[108]. Fortunately, it is possible to by-pass this result as 
explained in the next subsection. 



5.4.4 Prom Richardson-Gaudin to Korotkin-Nicolai equations 

In Ref.[109] bosonic and fermionic formalism for pairing models discussed in 
the previous subsection was developed. This formalism happens to be the most 
helpful for investigation of the gap problem. Indeed, define three operators 
= Y.m "am^im, = iA)^ ^ Uni AmArn ' ^hcy cau be used for construction 
of operators Ki = in/ ± jVli and K^^ ~ = (i^;~)^ such that they obey 
the following commutator algebra 

[Kl K+] = 5uK+, [K+,K-] - T'^SuKl (5.50) 

In this algebra as well as in the preceding expressions, the upper sign corresponds 
to bosons and the lower to fermions. In Ref. [79], we discussed such an algebra 
for the fermionic case only, e.g. see Eq.s (4.31) of [79]. These results can 
be extended now for the bosonic case. In fact, such an extension is already 
developed in Ref.[109]. Unlike [79], where we used sZ(2,C) Lie algebra, only 
its compact version, that is su(2), was used in [109] for representing fermions. 
For bosonic case the commutation relations, Eq.(5.50), are those for sm(1, 1) Lie 
algebra. Incidentally, in the paper by Korotkin and Nicolai, Ref. [31], exactly the 

^^The physical meaning of constants entering this equation is not important for us. It is 
given in Ref. [106].. 
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same Lie algebra was used. Furthermore, in the same paper it was argued that 
it is permissible to replace su{l, 1) by sl{2, R) Lie algebra while constructing the 
K-Z-type equations, e.g. read p. 428 of this reference. Since in [79] the sl{2, C) 
Lie algebra was used, that is a complexified version of sl{2. R), this allows us to 
use many results from our work. Thus, in this subsection we shall discuss only 
those results of [109] which are absent in our Ref.[79]. In particular, following 
this reference the set of Gaudin-like commuting Hamiltonians written in terms 
of operators Kf, and Kf is given by 

Hi = Kf + 2g{ ^ ^{K+K^ + KfK+) t Yu-KfKf,}. (5.51) 

Here Xw = Yw = [ei — ei/)"^ . For g ^ oo the first term can be ignored and the 
remainder can be used in the K-Z-type equations. The semiclassical treatment 
of these equations discussed in detail in [79] is resulting in the following set of 
Bethe (or R-G) ansatz equations 

n , m ^ 

to be compared with Eq.(5.47). Unlike Eq.(5.47), in the present case di = 
|(2i/i ±f];). The bosonic version of Eq. (5.52) corresponding to 5^(1, 1) Lie alge- 
bra coincides with Eq.(4.50) of Korotkin and Nicolai paper, Ref.[31], provided 
that the following identifications are made: di ^ si, 2ei jj. Unlike Ref.[31], 
where Eq.(5.52) was obtained by standard mathematical protocol, in this work 
it is obtained bascid on the underlying physics. Bccaiisc of this, it is appropriate 
to extend our physics-stype analysis by considering the case of finite g's. Then, 
Eq.(5.52) should be replaced by 

1=1 1=1 

In Ref.[31] the gap problem was discussed in detail for the fcrmionic case when 
the coupling constant g is negative (BCS pairing Hamiltonian) , e.g. see Eq.s 
(4.43)-(4.45) of Ref.[31]. In the present case we are dealing with the bosonic 
case for which the coupling constant is positive;. Hence the gap problem should 
be re analyzed. For this purpose, it is convenient to consider both positive and 
negative coupling constants in parallel for reasons which will become apparent 
upon reading. 

5.4.5 Emergence of the gap and the gap dilemma 

Eq.s(5.53) cannot be solved without some physical input. Initially, such an 
input was coming from nuclear physics (e.g. read [110-112] for general informa- 
tion on nuclear physics). Indeed, Richardson's papers [104,105] were written 
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having applications to nuclear physics in mind. Given this, the question arises 
about the place of the R-G model among other models describing nuclear spec- 
tra and nuclear properties. We need an answer to this question to finish proof 
of the gap's existence in QCD. 

Looking at the Gaudin-like Hamiltonian, Eq.(5.51), and comparing it with 
the Hamiltonian, Eq.(6), in Ref.[113p^ it is easy to notice that they are almost 
the same! Because of this, it becomes possible to transfer the methodology of 
Ref.[113] to the present case. Thus, it makes sense to recall briefly circumstances 
at which the gap emerges in nuclear physics. 

As is well known, the nuclei are made of protons and neutrons. One can 
talk about the number Af of nucleons, the number Z of protons and the number 
N of neutrons in a given nucleus. Nuclear and atomic properties happen to 
be interrelated. For instance, in analogy with atomic physics one can think of 
some effective nuclear potential in which nucleons can move "independently". 
This assumption leads to the shell model of nuclei. Use of Pauli principle guides 
fillings of shells the same way as it guides these fillings in atomic physics. This 
leads to emergence of magic numbers 2, 8, 20, 28, 50, 82 and 126 for either 
protons or neutrons for the totally filled shells. Accordingly, the most stable 
are the doubly magic (for both protons and neutrons) nuclei. It is of interest to 
know what kinds of excitations are possible in such shell models? The simplest 
of these is when some nucleon is moving from the closed shell to the empty shell 
thus forming a hole. When the number of nucleons increases, the question about 
the validity of the shell model emerges, again in analogy with atomic physics. As 
in atomic physics, one can think about the Hartree-Fock (H-F) and other many- 
body computational schemes, including that developed by Richardson- Sherman 
and Gaudin. For our purposes, it is sufficient to use only the Tamm-Dankoff 
(T-D) approximation to the H-F equations described, for example, in Ref.[112]. 
The essence of this approximation lies in restricting the particle-hole interac- 
tions to nucleons lying in the same shell. The T-D approximation is obtainable 
from the is R-G Eq.s(5.53) when the last term (effectively taking care of Pauli 
principle) in these equations is dropped. The T-D approximation was success- 
fully applied for description of the giant nuclear dipole resonance [110-112]. At 
the classical level the physics of this resonance was explained in the paper by 
Goldhaber and Teller [114]. The resonance is caused by two nuclear vibrational 
modes: one, when protons and neutrons move in the opposite directions and 
another- when they move in the same direction. Upon quantization of such 
classical model and taking into account the isotopic spin of nucleons, the trun- 
cated Eq.s(5.53) are obtained in which both signs for the coupling constant are 
allowed since the nucleon system is expected to be in two isospin states : T = 1 
and T — (@. Details of these calculations are given in Ref.[112], page 221. So- 
lutions of the T-D equations can be obtained graphically in complete analogy 
with that described in our work, Ref.[79]. These graphical solutions reflect the 
particle-hole duality built into the T-D approximation. Because of this duality, 

28published in 1961! 

^^This can be easily understood based on the fact that isospin for both particles and holes 
is equal to 1/2 [110-112]. 
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the magnitude of the gap in both cases should be the same. To demonstrate 
this, the seniority scheme described in [110-112] is helpful The seniority opera- 
tor was defined by Eq.(5.43). It determines the number of unpaired particles in 
the nuclear system. Since it commutes with the Hamiltonian, the many-body 
states can be classified with help of its eigenvalues vf. Suppose at first that all 
single particle energies £/ are the same (that is e/ = e) so that all seniority 
eigenvalues i^f are v. Let then Af be the total number of nucleons. Thus, the 
state for which i/ = contains only pairs, analogously, the state ly = 1 contains 
just one unpaired nucleon, v = 2 has 2 unpaired nucleons and J\f should be 
even and so on. So, states = 0,iy = 2,iy = 4, ... can exist only in even nuclei. 
For such nuclei the gap is nonzero. To see this, we follow Refs. [110-112] which 
we would like now to superimpose with the results of the Richardson-Sherman 
paper, Ref.[105]. Specifically, on page 231 of this reference one can find the 
following result for the ground state {v = 0) energy 

E^^o{N)=2Ne-gN{n-N+l) (5.54) 

where TV is the number of pairs. To connect this result with that in Refs. [110- 
112], let N = M 12 and consider the difference 

£^^q{N) = £^=0(^/2) -Me = -|aA(20 -M^2). (5.55) 

The obtained result coincides with Eq.(11.14) of Ref.[112] as required. To obtain 
states of seniority = 2n we use Eq.(3.2) of Ref.[105]. It reads 

E^^2n{N) = 2Ns - g {N - n) {n - N - n + 1), n = 0, N. (5.56) 

Repeating the same steps as in = case we obtain, 

g^{j^)=-l(j\f.i.)(^2n-Af-i^ + 2). (5.57) 
Finally, consider the difference 

£^{JV)-£^=oiJ^)=^v{2n-u + 2). (5.58) 

This result is in accord with Eq. (11.22) of Ref.[112]. Since the obtained dif- 
ference is A/'-independent it can be used both ways: a) for calculations in the 
thermodynamic limit A/" — > 00 and b) for making accurate calculations in the 
opposite limit of very small number of nucleons. In the simplest case we should 
consider only one shell and the first excited state of seniority 2 for this shell. 
Initially (the ground state) we have just one pair while finally (the first excited 
state) we have two independent particles occupying single particle levels. 

Looking at Eq.s(5.53) and letting there m = l(one pair) we recognize that 
the second sum in this set of equations disappears. Thus, by design, we are left 
with the T-D approximation. Using Eq.(5.58) for v — 2 we obtain the following 
value of the gap A : 

A = £2{J\f) - £o{Af) = gft. (5.59a) 
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Notice, that since is the degeneracy, there could be no more than Af = fl 
particles at the single particle level. Thus, in general we should have Af < Q. 
Because of the particle-hole duality, it is permissible to look also at the situation 
for which J\f > fJ.This is equivalent to changing the sign in front of the coupling 
constant. Repeating again all steps leads to the final result for the gap 



It is demonstrated in Ref.s [110-112] that in the limit A/" — >■ oo, when the contin- 
uum approximation (replacing summation by integration) can be used leading 
to a more familiar BCS-type equation for the gap, the result just obtained sur- 
vives. Indeed, in Ref.[103] the BCS-type result is obtained in the continuum 
approximation for the attractive Bose gas. In view of the results just obtained, 
it should be clear that such a result should hold for both attractive and repul- 
sive Bose gases. This conclusion is in accord with accurate recent Bethe ansatz 
calculations done in Ref.[115] for systems of finite size. Thus, we just arrived 
at the issue which we shall call tie gap dilemma. While the results obtained 
above strongly favor use of the repulsive Bose gas model (not linked with the 
F-S model), the results obtained in this subsection indicate that, after all, the 
F-S model ( linked with the attractive Bose gas model ) can also be used for de- 
scription of the ground and excited states of pure Y-M fields. The essence of 

THE DILEMMA LIES IN DECIDING WHICH OF THESE RESULTS SHOULD ACTUALLY 
BE USED. 

While the answer is provided in the next section, we are not yet done with 
the gap discussion. This is so because the seniority model is applicable only 
to the case when all single-particle levels have the same energy. This is too 
simplistic. We would like now to discuss more realistic case 

Before doing so, few comments are appropriate. In particular, with all suc- 
cesses of nuclear physics models, these models are much less convincing than 
those in atomic physics. Indeed, all nuclei are made of hadrons which are made 
of quarks and gluons. Thus the excitations in nuclei are in fact the excitations 
of quark-ghion plasma. This observation qualitatively explains why the R-G 
equations work well both in nuclear and particle physics. Some attempts to 
look at the processes in nuclear physics from the standpoint of hadron physics 
can bo foimd in Rcfs.[116,117]. 

Now we can return to the discussion of the T-D equations. Fortunately, de- 
tailed analytical study of these equations was recently made in Ref.[118]. The 
same authors extended these results to the case of two pairs in [119]. Since 
the results obtained in [119] are in qualitative agreement with those obtained 
in Ref.[118], we shall focus attention of our readers only on results of Ref.[118]. 
Thus, we need to find some kind of analytic solution of the following T-D equa- 
tion 



For different e-s normally it should have L eigenvalues {1 < fi < L). Since 
we are interested in finding the gap, the above equation is written for just one 



A = £2(AA)-fo(AA)=|g|f^- 



(5.59b) 




L 



(5.60) 
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nuclcon pair. Thus the seniority ^ = 0. It is of interest to check first what 
happens when ah e'^s coalesce. In such a case we obtain, 

fl 1 , , 

(5.61) 



2e- E g 

where f2 = fij and = e Vi = 1, L. Eq.(5.61) can be equivalently rewrit- 
ten as 

EQ = 2e- Qg. (5.62) 

This resuh for the ground state is in agreement with Eq.(5.54) for TV = 1. The 
first excited state is made of one broken pair so that the pairing disappears 
and the energy E^^2 = 2e. Prom here, the value of the gap is obtained as 
Ei,=2— E{) = gCl in agreement with Eq.(5.59). If now we make all energy 
levels different, then one can see that solutions to Eq.(5.60) are subdivided 
into those lying in between the single particle levels {trapped solutions) and 
those which lie outside these levels {collectivized solutions). For I9I sufficiently 
large the solution, Eq.(5.61), is the leading term (in the sense described below) 
representing the collectivized solution. Since the trapped solutions represent 
corrections to energies of single particle states, they do not contribute directly 
to the value of the gap. They do contribute to this value indirectly. Indeed, 
following Ref.[118] we rewrite Eq.(5.60) as 

X — ^ 1 V — ^ 1 

and expand the denominator of Eq.(5.63) in a power series. As result, the 
following expansion 

^^-^ = -1 - + ja^ + 0{a'') (5.64) 



is obtained in which e = j^J^^i^iy ^ — — /n^I^^il^i — e)"^ and 7 is 



related to the higher order moments ( details are in Ref.s[118,119]). Using these 
results, the gap is obtained in the same way as before. 

The quality of computations in Ref.[118] was tested for 3-dimensional har- 
monic oscillator (by adjusting dimensionality of this oscillator it can be thought 
of as "closed string model" representing both the shell model for atomic nu- 
cleus and the gluonic ring for the Y-M fields) for which Si = {i + 3/2) (in 
the system of units in which hb: = \) and fij = (i -|- l)(i -I- 2)/2. For this 3- 
dimensional oscillator corrections to the collectivized energy, Eq.(5.64), become 
negligible already for \g\ > 0.2, provided that L > 8. Obtained results allow us 
to close this section at this point. These results are of no help in solving the 
gap dilemma though. This task is accomplished in the next section. 
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6 Resolution of the gap dilemma 



6.1 Motivation 

In the previous section we provided evidence linking the gap problem for Y-M 
fields with the problem about the excitation spectrum of the repulsive Bose gas. 
The gap equation, Eq.(5.59), is also used in nuclear physics where it is known 
to produce the same value for the gap for both signs of the coupling constant g. 
Since both options are realizable in Nature in the case of nuclear physics, the 
question arises about such possibility in the present case. In the case of nuclear 
physics experimental realization (giant nuclear dipole resonance) of both options 
for the coupling constant is experimentally testable. Thus, in the present case 
we have to find some alternative physical evidence. If, indeed, such evidence 
could be found, this would allow us to bring back into play the well studied F-S 
model which microscopically is essentially equivalent to the XXX Id Heisenberg 
ferromagnet as results of Appendix B and subsections 3.5 and 5.2 indicate. The 
next subsection supplies us with the alternative physical evidence. 

6.2 Some facts about harmonic maps and their uses in 
general relativity 

Suppose we are interested in a map from m— dimensional Riemannian manifold 
M with coordinates a;" and metric 7^j(x) to n-dimensional Riemannian man- 
ifold A/" with coordinates and metric GAsif)- A map — > A/" is called 
harmonic if ip"^ (x") satisfies the Euler-Lagrange (E-L) equations originating 
from minimization of the following Lagrangian 

c ^ ^GABiv)r\^)vi^^b (6.1) 

in which 7 = det(7^j). Since such defined Lagrangian is part of the La- 
grangian given by Eq.(3.6), the E-L equations for Eq.(6.1), in fact, coincide 
with Eq.s(3.10). In the most general form they can be written as [38 

V'i'a^+r^c^^^'" =0. (6.2) 

In such a form we can look at transformations ip^ = ip^ {^^) keeping C form- 
invariant. To find such transformations, following Neugebauer and Kramer [38], 
we introduce the auxiliary Riemannian space defined by the metric 

dS^ = GAB{^)d'p'^dif^. (6.3) 

Use of the above metric allows us to investigate the invariance of C with help of 
standard methods of Riemannian geometry. In the present case, this means that 
one should study Killing's equations in spaces with metric Gab- Specifically, let 
us consider the Lagrangian for source-free Einstein-Maxwell fields admitting at 

•^"We use the 1st edition of Ref. [38] for writing this equation. This means that we have to 
define T^^ as T^^ = Ig"-^! gf^Gw + ^G,a - T^Gt.}. 
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least one non-null Killing vector ^. To design such a Lagrangian we begin with 
the Ernst equation, Eq.(2.4), for pure gravity and replace the Ernst potential 
e = —F + by two complex potentials £ and Then, by symmetry, the 
equations for stationary Einstein-Maxwell fields can be written as follows [38] 

F£:2 + 7°'f ,a(f 6 + 2$,bl') = 0, + -f^''^A£,b + = 0. (6.4) 

These equations are obtained by minimization of the Lagrangian 

C = Vl[Rab + 2F~^r''^^a^^b + \F-^l''\£,a + 2l>$,,)(f + 2$$^b)], (6.5) 

i.e. from equations = 0, |§ =0 and |^ = 0. Taking these results into 
account, the auxiliary metric, Eq.(6.3), can now be written as 

dS^ = 2F-^d<^d^ + i^-^ \d£ + 2$d$|^ . (6.6) 

The analysis done by Neugebauer and Kramer [38] shows that there are eight 
independent Killing vectors leading to the following finite transformations : 

£' = aa£, (f>' = a<i>; 

£' =£ + ib, $' = $; 

£' = £{l +_ic£)-\ $' = (1 + ic£)-i; (6.7) 

= £ - 2/3$ - /3/3, $' = $ -f ;3; 

£' = £{1 - 27$ - 775)-\ $' = ($ + 7£)(1 - 27$ - j^£)-\ 

Complex parameters a,P,^ as well as real parameters b and c are connected 
with these eight symmetries. Evidently, solutions f are also solutions of 
Eq.s(6.4), provided that 7"^ stays the same. Therefore if, say, we choose some 
vacuum solution as a " seed" , we would obtain, say, the electrovacuum solution in 
accord with Appendix A. Incidentally, the electrovacuum solutions obtained by 
Bonnor (Appendix A) cannot be obtained with help of transformations given by 
Eq.s(6.7). They are considered separately below. These observations allow us to 
reduce the Lagrangian C to the absolute minimum without loss of information. 
In 1973 Kinnersley [38] found that the group of symmetry transformations for 
the Einstein-Maxwell equations with non null Killing vector is the group SU(2,1) 
which has eight independent generators. In view of the above mentioned reduc- 
tion of C it is sufficient to replace the metric in Eq.(6.6) by a collection of much 
simpler metric related to each other by transformations Eq.(6.7). All the possi- 
bilities are described in the Table 34.1 of Rcf.[38]. For our needs we focus only 
on three of these (much simpler/reduced) metric listed in this table. These are 

dS^ = -m-.^l^, (6.8) 
- (i_ci>$)2' ^^-^^ 



^^Recall, that —F = V according to notations introduced in connection with Eq.(2.4). 



45 



and 



(l + $$)2 

The first and the second of these metric, corrc^spond to the vacuum state, 
respectively, with $ = and £ = —1, of pure gravity associated with the 
subgroup SU(1,1) of SU(2,1). The third metric, Eq.(6.10), corresponds to a 
subgroup SU(2). It is related to the electrostatic fields {£ = 1) such that the 
space-time becomes asymptotically flat for >O.It is important that the metric, 
Eq.(6.10), is related to the vacuum metric, Eq.s(6.8),(6.9), via transformations 
either listed in Eq.(6.7) or related to these transformations. In particular, the 
related transformations can be obtained as follows. Using Ref. [38] , it is conve- 
nient to make the parameters b and c in Eq.s(6.7) complex and to consider all 
eight complex parameters as independent of their complex conjugates. Under 
such conditions the metric given by Eq.(6.10) is related to that given by Eq.(6.8) 
by the simplest complex transformation: = and = . These transfor- 
mations indicate that, starting with real vacuum solution for pure gravity as a 
seed, the above transformations are capable of reproducing some electrovacuum 
solutions. Additional details are discussed below. 

These results can be interpreted as follows. While the Ernst functional, 
Eq.(3.18), is representing pure axially symmetric gravity, the F-S-type func- 
tional, Eq.(3.19), should describe some special case of electrovacuum (Maxwell- 
Einstein) gravity. In view of results of Appendix C, it is possible to use these 
transformations in reverse (see below), that is to obtain the results for pure 
gravity from those for electrovacuum. This peculiar " duality" property of grav- 
itational flelds provides physically motivated resolution of the gap dilemma and, 
in addition, it allows us to obtain many new results. 

6.3 Resolution of the gap dilemma and SU(3)xSU(2)xU(l) 
symmetry of the Standard Model 

The original F-S-type model thus far is limited only to SU(2) gauge theory. 
SU(2) gauge theory is known to be used for description of electroweak interac- 
tions where, in fact, one has to use the gauge group SU(2)xSU(l) [19]. The 
hadron physics (that is QCD) requires us to use the gauge group SU(3). This 
is caused by the fact that quark model of hadrons uses flavors (e.g. u,d,s,c,b, t) 
labeling quarks of different masses. Each of these quarks can be in three differ- 
ent colors (r,g,b) standing for "red", "green" and "blue". Presence of different 
colors leads to fractional charges for quarks. Far from the target the scattering 
products are always colorless. The gauge group SU(3) is used for description of 
these colors. Although theoretically the number of colors can be greater than 
three, this number is strictly three experimentally [19]. The results of this work 
allow us to reproduce this number of colors. For this purpose we have to be 
able to provide the answer to the following fundamental question: 



(6.10) 
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Can equivalence between gravity and Y-M fields (for SU(2) 

GAUGE group) DISCOVERED BY LOUIS WiTTEN BE EXTENDED TO THE GROUP 
SU(3)? 

Very fortunately, this can be done! For the; sake of space, we shall be brief 
whenever details can be found in literature, e.g. see Refs. [120-122]. 

To proceed, first, we have to go back to Eq.s(2.14),(2.15) and to modify 
these equations in such a way that instead of the Ernst Eq.(2.4) for the vacuum 
(gravity) field we should be able to obtain Eq.s (6.4) for electrovacuum. In 
the limit $ = the obtained set of equations should be reducible to Eq.(2.4). 
As it was noticed by Giirses and Xanthopoulos [120], in general, this task can- 
not be accomplished. Indeed, these authors demonstrated that the self-duality 
Eq.s(2.14) for SU(2) and for SU(3) Lie groups look exactly the same for axi- 
ally symmetric fields. Nevertheless, in the last case, upon explicit computation 
instead of the vacuum Ernst Eq.(2.4) one gets an electovacuum equations (e.g. 
see Eq.s(6.4)) which, following Ernst [43], can be explicitly written as 

(Re£ + |$|^) V^S = {VS + 2$V$) • VS, (6.11a) 
(Re£ + = (V£ + 2$V$) • V$. (6.11b) 



M = f- 



(6.12) 



These equations are obtained if, instead of the matrix M given by Eq.(2.15), 
one uses 

1 ^/2$ -1(5-5-2$$)" 

V2$ -1(5-1-5-2$$) -i\/2$5 

j(5-5-2$$) iv^5$ 5 5 

in which, instead of the one complex potential e = —F + used for solution of 
the vacuum Ernst Eq.(2.4), two complex potentials 5 and $ are being used. In 
this expression the overbars denote the complex conjugation and / = — ^ (e + 
e + 2$$). Since the Einstein-Maxwell Eq.s(6.4) (or (6.11)) are invariant with 
respect to transformations given by Eq.s(6.7), there should be a matrix A with 
constant coefficients such that the M' = AMA^ will have primed potentials 5 
and $ taken from those listed in the set Eq.(6.7). Authors of [120] found explicit 
form of such A-matrices. However, when instead of matrix M we substitute the 
matrix M' into self-duality Eq.s(2.14), the combination M'~^dM' looses this 
information. As result, we are left with the following situation: while on the 
gravity side the matrix M' = AMA^ does allow us to obtain new and physically 
meaningful solutions from the old ones, on the Y-M side all this information 
is lost. Thus, the one-to-one correspondence discovered by L.Witten for SU(2) 
is apparently lost for SU(3). Very fortunately, this happens only apparently! 
This is so because the Neugebauer- Kramer (N-K) transformations described 
by Eq.s(6.7) do not exhaust all possible transformations which can be applied 
to the matrix M, Eq.(6.12). Among those which are not accounted by N- 
K transformations are those by Bonnor [38,123] whose work is mentioned in 
Appendix A. These are given by 

5 = ee;$ = ^(e-e) = iw, (6.13) 
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where e = —F + zw is solution of the Ernst Eq.(2.4). In view of the results of 
Appendix A one can be sure that the potentials £ and $ satisfy Eq.s(6.11). 
This means that one can use these (Bonnor's) potentials in the matrix M to 
reproduce Eq.s(6.11). This time, there is onc-to one correspondence between the 
self-duality Y-M and the Einstein-Maxwell equations. Even though this is true, 
the question immediately arises about relevance of such solutions to the solution 
of the gap problem discussed in Section 5. In Section 5 the Ernst Eq.(2.4) was 
used essentially for this purpose while Eq.s(6.11) are seemingly different from 
Eq.(2.4). Again, fortunately, the difference is only apparent. 

From the definition of Bonnor transformations, Eq.(6.13), it follows that 
the potential £ is real. Also, from the same definition it follows that |$| = 
w^. Introduce now new potential Z = £ -{-uj^. For it, we obtain 

VZ ^ V{£ + Lu^) = V£+2ujVlo ^V£ + 2$V$. (6.14) 

Using this result, Eq.s(6.11) can be rewritten as follows 

(ZV^-VZ • V) ^ ^ = 0. (6.15) 

Furthermore, consider the related equation 

{ZV^-VZ • V) = 0. (6.16) 

Evidently, if it can be solved, then equation (ZV^— VZ- V)w = can be 
solved as well. This being the case, the system of Eq.s(6.15) will be solved if 
the Ernst-type vacuum equation 

ZV'^=VZVZ (6.17) 

of the same type as Eq.(2.4) is solved. The obtained result is opposite to that 
derived by Bonnor, described in Appendix A (see also works Hauser and Ernst 
[124] and by Ivanov [125]). This means that, at least in some cases (having 
physical significance) the self-dual Y-M fields for both SU(2) and SU(3) gauge 
groups are obtainable as solutions of the Ernst Eq.(2.4). This means that all 
results of Section 5 obtained for SU(2) go through for the gauge group SU(3). 

With these results at our disposal we would like to discuss their applications 
to the Standard Model [19, 126]. From Ref. [120] it is known that the matrix 
M e SU(3) has subgroups which belong to SU(2). In particular, one of such 
subgroups is obtained if we let $ = in Eq.(6.12). Then, in view of Eq.(6.17), 
it is permissible to replace f by e of Eq.(2.4). Thus, the obtained matrix M is 
decomposable as M = Mi + M2, where the matrix Mi is given by 



Ml = /-I 



1 CO 

w ee 



(6.18) 
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in agreement with the matrix M defined by Eq.(2.15) since in this case / 
— + e) = F. At the same time, the matrix is given by 



Mo 




1 




(6.19) 



Using elementary operations with matrices we can represent matrix M in the 
form 

"0 1" 

M= a 6 (6.20) 
6 c 

where a = 1/F, b = lu/F and c = (F^ + lu'^)/F. Such a form of the matrix 
M is typical for the semidirect product of groups (when group elements are 
represented by matrices). In general case one should replace M by 



M = 



1 
a b ai 
b c a2 



Since the 2x2 submatrix belongs to SU(2) (because its determinant is 1) nor- 
mally describing a rotation in 3d space (in view of SU(2)<=^SO(3) correspon- 
dence), the parameters ai and 02 are responsible for translation. In this, more 
general case, the matrix AI describes the Galilean transformations, that is a 
combination of translations and rotations. If the translational motion is one 
dimensional it can be compactified to a circle in which case we obtain the cen- 
tralizer of SU(3) as SU(2) xU(l). At the level of Lie algebra su(3) this result was 
obtained in Ref.[127], pages 232 and 267. Its physical interpretation discussed 
in this reference is essentially the same as ours. The obtained centralizer is 
the symmetry group of the Wcinberg-Salam model (part of the standard model 
describing electroweak interactions). 

All these arguments were meant only to demonstrate that the F-S-type 
model, Eq.(3.19), should be used for description of electroweak interactions. 
For description of strong interactions, in accord with Ref.[120], we claim that 
the matrix M given by Eq.(6.12) in which £ and $ are taken from Bonnor's 
Eq.s(6.13) is intrinsically of SU(3) type. That is, it cannot be obtained from 
the matrix M (in which = 0) by applications of the N-K transformations, 
i.e. there are no transformations of the type M'($) = AM{^ = 0)A^ . There- 
fore, this type of SU(3) matrix should be associated with QCD part of the SM. 
Hence we have to use the Ernst functional, Eq.(3.18), instead of the F-S-type, 
Eq.(3.19). These results provide resolution of the gap dilemma. Evidently, this 
resolution is equivalent to the statement that the symmetry group of the SM is 
SU(3)xSU(2)xU(l). This result should be taken into account in designing all 
possible grand unified theories (GUT). In the next subsection we shall discuss 
the rigidity of this result. 
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6.3.1 Remarkable rigidity of symmetries of the Standard Model and 

the extended Ricci flow 

In addition to Bonnor's transformations there are many other transformations 
from vacuum to electrovacuum. In particular, in Appendix A we mentioned 
transformations discovered by Hcrlt. By looking at Eq.s(A.5)-(A.7) describing 
these transformations and comparing them with those by Bonnor, Eq.(6.13), 
it is an easy exercise to check that all arguments leading from Eq.s(6.11) to 
(6.17) go through unchanged. By using superposition of N-K transformations 
and those either by Bonnor or by Herlt it is possible to generate a countable 
infinity of vacuum-to electrovacuum transformations such that they could be 
brought back to the vacuum Ernst solution, Eq.(6.17), using results of previous 
subsection. This property of Einstein and Einstein -Maxwell equations we shall 
call "rigidity". In view of results of previous subsection, this rigidity explains 
the remarkable empirical rigidity of symmetries of the SM. Indeed, suppose that 
the color subgroup SU(3) can be replaced by SU(N), N>3. In such a case it 
is appropriate again to pose a question : Can self-dual Y-M fields-gravity cor- 
respondence discovered by L.Wittcn for SU(2) be extended for SU(N), N>3? 
In Ref.[128] Giirses demonstrated that, indeed, this is possible but under non- 
physical conditions. Indeed, this correspondence requires for SU(n-|-l) self-dual 
Y-M fields to be in correspondence with the set of n-1 Einstein-Maxwell fields. 
Since n = 1 and n = 2 cases have been already described, we need only to 
worry about n > 2. In such a case we shall have many-to-one correspondence 
between the replicas of electrovacuum and vacuum Einstein fields which, while 
permissible mathematically, is not permissible physically since the Bonnor-type 
transformations require one-to-one correspondence between the vacuum and 
electrovacuum fields. Herrcra-Aguillar and Kechkin, Ref.[129], found a way of 
transforming the compactified fields of heterotic string (e.g. see Eq.(3.12)) into 
Einstein-multi-Maxwell fields of exactly the same type as discussed in the paper 
by Giirses [128]. While in the paper by Giirses these replicas of Maxwell's fields 
needed to be postulated, in [129] their stringy origin was found explicitly. From 
here, it follows that results obtained in this subsection make the minimal func- 
tional, Eq.(3.8), and the associated with it Pcrclman-like functional, Eq.(3.13), 
universal. The universality of the associated with it Ricci flow, Eq.s (3.14), 
has physical significance to be discussed below. 

7 Discussion 

7.1 Connections with loop quantum gravity 

A large portion of this paper was spent on justification, extension and exploita- 
tion of the remarkable correspondence between gravity and self-dual Y-M fields 
noticed by Louis Witten. Such correspondence is achievable only non per- 
turbatively. In a different form it was emphasized in the paper by Mason and 
Newman [130] inspired by work by Ashtekar, Jacobson and Smolin [131]. It is 
not too difficult to notice that, in fact, papers [130,131] are compatible with 
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Witten's result since reobtaining of Nahm's equations in the context of grav- 
ity is the main result of Ref.[131]. In this context the Nahm equations are 
just equations for moving triad on some 3-manifold. Since the connection of 
Nam's equations with monopoles can be found in Ref.[68] and with instantons 
in Ref.[132] the link with Witten's results can be established, in principle. Since 
the authors of [131] are the main proponents of loop quantum gravity (LQG) 
such refinements might be helpful for developments in the field of LQG. We 
shall continue our discussion of LQG in the next subsection. 



7.2 Topology changing processes, the extended Ricci flow 
and the Higgs boson 

According to the existing opinion the SM does not account for effects of grav- 
ity. At the same time, in the Introduction we mentioned that in recent works 



by Smolin and collaborators [32-34 
certain quantum gravity theorie^ 



it was shown that "topological features of 
can be interpreted as particles, matching 
known fermions and bosons of the first generation in the Standard Model" . 
Similar results were also independently obtained in works by Finkelstein, e.g. 
see Ref.[133] and references therein. In particular, Finkelstein recognized that 
ah quantum numbers describing basic building blocks(=particles) of the SM can 
be neatly organized with help of numbers used for description of knots. More 
precisely, with projections of these knots onto some plane. It happens, that for 
description of all particles of the electroweak portion of the SM the numbers 
describing trefoil knot are sufficient. The task of topological/knotty description 
of the entire SM was accomplished to some extent in Ref.[33]. This reference 
as well as Ref.s[32-34] in addition are capable of describing particle dynam- 
ics/transformations. All these works share one common feature: calculations 
do not require Higgs boson. This fact is consistent with results discussed in 
subsection 4.3.1. 

The question arises: Is this feature a serious deficiency of these topological 
methods or are these methods so superior to other, that the Higgs boson should 
be looked upon as an artifact of the previously existing perturbative methods 
used in SM calculations? To answer this self-imposed question requires several 
steps. 

First, we recall that according to the existing opinion the SM does not 
account for effects of gravity. In such a case all the above results should have 
nothing in common with the SM which is not true. 

Second, the results obtained in this paper indicate that knots/links/braids 
mentioned above have not only virtual (combinatorial/topological) but also 
differential-geometric description (Appendix B). Because of this, topological 
description should be looked upon as complementary to that obtainable with 
help of the F-S-type models. 

Third, it is known that knot/link- describing Faddeev model can be con- 
verted into Skyrme model [134]. It is also known that the Skyrme-type models 



32That is LQG. 
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do not account for quarks explicitly, Ref.[68], page 349. This is not a serious 
drawback as we shall explain momentarily. 

Fourth, much more important for us is the fact that the Skyrme model can 
be used both in nuclear [135] and high energy [136] physics where it is used for 
description of both QCD (nicely describing the entire known hadron spectra) 
and electroweak interactions. 

To account for quarks one has to go back to the Faddecv-type models capa- 
ble of describing knots/links and to make a connection between these physical 
knots/links and topological/combinatorial knots/ links discussed in Refs[32- 
34,133]. This is still insufficient! It is insufficient because Floer's Eq.(4.7) con- 
nects different vacua each is being described by the zero curvature condition 
Eq.(4.13). It is always possible to look at such a condition as describing some 
knot /link differential geometrically. With each knot, say in S^, some 3-manifold 
is associated. Furthermore such a manifold should be hyperbolic (subsection 
3.6), that is either associated with hyperbolic-type knot/links [20,137] in 
or with knots/links "living" in hyperboloid embedded in the Minkowski space- 
time. Such a restriction is absent in Ref.s[32-34,133]. At the same time the Y-M 
functional, Eq.(4.12), is defined for a particular 3-manifold whose construction 
is quite sophisticated. Eq.(4.7) describes processes of topology change by con- 
necting difi'erent vacua. Such changes formally are not compatible with the fact 
that we are dealing with one and the same 3-manifold Mx [0,1]. Prom the math- 
ematical standpoint [11] no harm is made if one considers just this 3-manifold, 
e.g. read Ref.[ll], page 22, bottom. Since particle dynamics is encoded in 
dynamics of transformations between knots/links, it causes us to consider tran- 
sitions between different 3-manifolds. These 3-manifolds should be carefully 
glued together as described in Ref.[ll]. In this picture particle dynamics involv- 
ing particle scattering/transformation is synonymous with processes involving 
topology change. These are carried out naturally by instantons. Such processes 
can be equivalently and more physically described in terms of the properties of 
the (extended) Ricci flow (subsection 3.4) following ideas of Perelman's proof of 
the Poincare' conjecture. Indeed, experimentally there is only finite number of 
stable particles. Without an exception, the end products of all scattering pro- 
cesses involve only stable particles. This observation matches perfectly with the 
irreversibility of Ricci flow processes involving changes in topology: from more 
complex-to less complex 3-manifolds. Such Ricci flow model upon development 
could provide mathematical justification to otherwise rather vague statements 
by Finkelstein that " more complicated knots ( particles) can therefore dynami- 
cally decay to trefoils (stable particles)", Ref.[133], page 10, bottom. 

7.3 Elementairy pcirticles as black holes 

In the paper [138] by Reina and Treves and also in [139] by Ernst it was found 
that for asymptotically flat Einstein-Maxwell fields generated from the vacuum 
fields by means of transformations of the type described above, in Section 6, 
the gyromagnetic factor g = 2. For the sake of space, we refer our readers to a 
recent review by Pfister and King [140] for definitions of g and many historical 
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facts and developments. In [140] it was noticed that such value of g is typical 
for most of stable particles of the SM. In view of the quantum gravity- Y-M 
correspondence promoted in this paper, the interpretation of elementary par- 
ticles as black holes makes sense, especially in view of the following excerpt 
from Ref.[38], page 526, "There is one-to-one correspondence between station- 
ary vacuum fields with sources characterized by masses and angular momenta 
and stationary Einstein-Maxwell fields with purely electromagnetic sources, i.e. 
charges and currents." 

Appendix A 

Peculiar interrelationship between gravitational, electromagnetic 
and other fields 

Unification of gravity and electromagnetism was initiated by Nordstrom in 
1913- before general relativity was formulated by Einstein. Almost immediately 
after Einstein's formulation, Kaluza, in 1921, and Klein, in 1926, proposed uni- 
fication of electromagnetism and gravity by embedding Einstein's 4-dimensional 
theory into 5 dimensional space in which 5th dimension is a circle. These results 
and their generalizations (up to 1987) can be found in the collection of papers 
compiled by Applequist, Chodos and Preund [141]. Regrettably, this collection 
does not contain alternative theories of unification. Since such alternative theo- 
ries are much less known/popular to/with string and gravity theoreticians, here 
we provide a brief representative sketch of these alternative theories. 

The 1st unified Einstein-Maxwell theory in 4-dirncnssional space-time was 
proposed and solved by Rainich in 1925. It was discussed in great detail by 
Misner and Wheeler [142]. After Rainich there appeared many other works on 
exact solutions of Einstein-Maxwell fields [38]. The most striking outcome of 
these, more recent, works is the fact that multitude of exact solutions of the 
combined Einstein-Maxwell equations can be obtained from solutions of the 
vacuum Einstein equations. 

In 1961 Bonnor [123]obtained the following remarkable result (e.g. read his 
Theorem 1). Suppose solutions of the vacuum Einstein equations are known. 
Using these solutions, it is possible to obtain a certain class of solutions of 
Einstein-Maxwell equations. 

In Section 6 we obtained the reverse result: Einstein's solutions for pure 
gravity were obtained from solutions of the Einstein-Maxwell cqiiations. With- 
out doing extra work, the electrovacuum solution obtained by Bonnor can be 
converted into that describing propagation of the combined cylindrical gravita- 
tional and electromagnetic waves. With some additional efforts one can use the 
obtained results as an input for results describing the combined gravitational, 
electromagnetic and neutrino wave propagation [143-144]. 

The results by Bonnor comprise only a small portion of results connecting 
static gravity fields with electromagnetic and neutrino fields. The next example 
belongs to Herlt [38,145]. It provides a flavor of how this could be achieved. 
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We begin with Eq.(2.5). When written expUcitly, this equation reads 

(^l + -^dp + dl^u = o. (A.i) 

This type of sohition is the result of use of the matrix M, Eq.(2.15), in Eq.(2.14b). 
Nakamura [146] demonstrated that there is another matrix Q given by 

^ = ( /i pu^^-p^r^ ) (^-^^ 

and the associated with it analog of Eq.(2.14b) 

dp{pdpQ ■ g-i) + d,{pd,Q ■ Q-i) = (A.3) 
leading to the equation analogous to Eq.(A.l), that is 

[dl-U, + di^u = Q. (A.4) 

Nakamura demonstrated that the solution u is obtainable from solution of 
Eq.(A.l). and vice versa. Thus, instead of the Ernst Eq.(2.4) we can use 
Eq.(A.4). This fact plays crucial role in Hertl's work. In it, he uses Eq.(A.4) 
to obtain u in Eq.(A.l) as follows 

exp(2u) = (u^i +G)^ (A.5) 

with G given by 

G = u,p[pi'>J%p + «fJ - uu,p]~\ (A.6) 
These results allow him to introduce a potential x via 

X = u-'-G. (A.7) 

Using the original work of Ernst [43] as well as Ref. [38] , we find that solution 
of the static axially symmetric coupled Einstein-Maxwell equations is given in 
terms of complex potentials e and In particular, in purely electrostatic case 
one has e = e =e^" — x and ^ = ^ = x while the magnetostatic case is obtained 
from the electrostatic by requiring = $ = i/; and e = e =e^" — . In this 
case is just relabeled x- Ref. [38] contains many other examples of the cou- 
pled Einstein-Maxwell equations obtained from the vacuum solutions of Einstein 
equations. 

The above results should be looked upon from the standpoint of fundamental 
problem of the energy-momentum conservation in general relativity requiring in- 
troduction (in the simplest case) of the Landau-Lifshitz (L-L) energy-momentum 
pseudotensor. The description of more complicated pseudotensors (incorporat- 
ing that by L-L) can be found in the monograph by Ortin [147]. To this one 
should add the problem about the positivity of mass in general relativity. The 
difficulties with these concepts stem from the very basic observation, lying at the 
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heart of general relativity, that at any given point of space-time gravity field 
can be eliminated by moving in the appropriately chosen accelerating frame 
(the equivalence principle) . This fact leaves unexplained the origin of the tidal 
forces requiring observation of motion of at least two test particles separated 
by some nonzero distance. The explanation of this phenomenon within general 
relativity framework is nontrivial.lt can be found in [148]. In turn, it leads to 
speculations about the limiting procedure leading to elimination of gravity at a 
given poiniF^. Apparently, this problem is still not solved rigorously [147]. An 
outstandig collection of rigorous results on general relativity can be found in 
the recent monograph by Choquet-Bruhat [149] while [150] discusses peculiar 
relationship between the Newtonian and Einsteinian gravities at the scale of 
our Solar system. 

Conversely, one can think of other fields at the point/domain where gravity 
is absent as subtle manifestations of gravity. Interestingly enough, such an idea 
was originally put forward by Rainich already in 1925 ! Recent status of these 
ideas is given in paper by Ivanov [125]. From such a standpoint, the functional 
given by Eq.(3.13) (that is the Perelman-like entropy functional) is sufficient for 
description of all fields with integer spin. With minor modifications (e.g. in- 
volving either the Newman-Penrose formalism [143,144] or supersymmeric for- 
malism used in calculation of Seiberg-Wittcn invariants [66]), it can be used for 
description of all known fields in nature. 

Appendix B 

Some facts about integrable dynamics of knotted vortex filaments 

B.l Connection with the Landau-Lifshitz equation 

Following Ref. [85] , we discuss motion of a vortex filament in the incompress- 
ible fluid. Some historical facts relating this problem to string theory are given 
in our recent work, Ref. [84]. Let u be a velocity field in the fluid such that 
divu = 0. Therefore, we can write u = V x A . Next, we define the vorticity 
w = V X u so that eventually, 



This expression can be simplifled by assuming that there is a Une vortex which 
is modelled by a tube with a cross-sectional area dA and such that the vorticity 
w is everywhere tangent to the line vortex and has a constant magnitude w. 
Let then F = J wdA so that 



with df being an infinitesimal line segment along the vortex. Such a model 
of a vortex resembles very much model used for description of dynamics of 
ring polymers [84]. Because of this, it is convenient to make the following 

'^''The abundance of available energy- momentum pseudotensors is result of these specula- 
tions. 




(B.la) 




(B.lb) 
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identification : u{'y{s,t)) = ^{s,t), with s being a position along the vortex 



at 

contour and i-time. This allows us to write 



^(s' t)- ^ i (7(^7) -7GM)) dj^ 
dt^ ' ' AttJ \\j{s',t)-j{s,t)f ds 



and to make a Taylor series expansion in order to rewrite Eq.(Blc) as 
97 r dj 9^7 f ds 



dt in'ds' ^ 5s'2 



+ ...]. (B.ld) 



In this expression only the leading order result is written explicitly. By intro- 
ducing a cut off £ such that \s — s'\> e and by rescaling time: t — > ^fln(£~^) 
one finally arrives at the basic vortex filament equation 

^7 _ 97 

Introduce now the Serret-Frenct frame made of vectors B,T and N so that 
B = T X N, kN T = 1^, where k is a curvature of 7. Then, Eq.(B.2) can 
be equivalently rewritten as 

|^ = .B (B.3) 

or, as 

dT 

— = T X T,,. (B.4) 

In the last equation the replacement s ^ x was made so that the obtained 
equation coincides with the Landau-Lifshitz (L-L) equation describing dynamics 

of Id Heisenberg ferromagnets [86]. 

B.2 Hashimoto map and the Gross-Pitaevskii equation 

Hashimoto [85] found ingenious way to transform the L-L equation into the 
nonlinear Scrodinger equation (NLSE) which is also widely known in condensed 
matter physics literature as the Gross-Pitaevskii (G-P) equation [96]. Because 
of its is uses in nonlinear optics and in condensed matter physics for descrip- 
tion of the Bose-Einstein condensation (BEC) theory of this equation is well 
developed. Some facts from this theory are discussed in the main text. Here we 
provide a sketch of how Hashimoto arrived at his result. 

Let T,U and V be another triad such that 

XX XX 

U = cos( j Tds)N - sin( j rds)B, V = sin( j rds)N + cos( J rds)B (B.5) 

in which r is the torsion of the curve. Introduce new curvatures /ti and K2 in 
such a way that 

X X 

Ki = KCOs( J rds) and K2 = ksm{ J rds), 
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then, it can be shown that 

^ = -K2U + KiV (B.6a) 
at 

and 

-J-=KiV + K2V. (B.6b) 

ox/ 

Using these equations and taking into account that UtV = — UVt after some 
algebra one obtains the following equation 

^V't + ^xx + [^l^l'-^W]V' = (B.7) 

in which ij) = ki+ iK2 and A{t) is some arbitrary x-independent function. By 

t 

replacing with V'exp(— i J dt'A{t')) in this equation we arrive at the canonical 
form of the NLSE which is also known as focussing cubic NLSE. 

+ ^ 1^1^ = (B.8) 

It can be shown that its solution allows us to restore the shape of the curve/filament 
7(5, t). The G-P equation can be identified with Eq.(B.7) if we make A{t) time- 
independent. In its canonical form it is written as (in the system of units in 
which ft = 1, m = 1/2) [86] 

iVt = -Vxx + 2K(|Vl'-c')V' = 0. (B.9) 

In general, the sign of the coupling constant k can be both positive and negative. 
In view of Eq.(B.8), when motion of the vortex filament takes place in Euclidean 
space, the sign of k is negative. This is important if one is interested in dynamic 
of knotted vortex filaments [85] . For purposes of this work it is also of interest to 
study motion of vortex filaments in the Minkowski and related (hyperbolic, de 
Sitter ) spaces. This should be done with some caution since the transition from 
Eq.(B.la) to (B.2) is specific for Euclidean space. Thus, study can be made at 
the level of Eq.s (B.3) and (B.4). Fortunately, such study was performed quite 
recently [94,95]. The summary of results obtained in these papers can be made 
with help of the following definitions. Introduce a vector n = {ni,n2,n3} so 
that the unit sphere S"^ is defined by 

S"^ :nl+nl+nl = 1. (B.IO) 

Respectively, the de Sitter space S^'^ (or unit pseudo sphere in R^'^ ) is defined 
by 

S^'^ : nj + nj ~ nl = 1, (B.ll) 
while the hyperbolic space H^(or hypcrboloid embedded in R^'^) is defined by 

n"^ :nl + nl-nl = -l,n3> 0. (B.12) 
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Using these definitions, it was proven in [94,95] that: a) for both de Sitter and 
spaces there are analogs of the L-L equation ( e.g. those discussed in the 
main text, in subsection 5.2.); b) the Hasimoto map can be extended for these 
spaces so that the respective L-L equations are transformed into the same NLSE 
(or G-P) equation in which k is positive. 
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